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ABSTR ACT 


Radar  is  considered  as  a  means  of  producing  a  map  in  two 
coordinates  (x  and  y)  of  a  dense  distribution  of  scatterers. 

It  is  assumed  that  a  complex  un it-scatterer  return  h(t;x,y) 
and  a  complex  scatterer-dens  ity  function  \Jr(x,y)  can  be  de¬ 
fined  such  that:  (a)  hCt  jx,y)  depends  only  on  the  radar; 

(b)  >Jr(x  ,y )  depends  only  on  the  scatterer  distribution;  (c) 
the  complex  video  representation  of  the  echo  signal  is 

r00 

s(t)  =  /  i)r(x,y)  h(tjx,y)  dxdy;  and  (d)  an  x-y  display  of 
*'—CO 

i|r(x,y)  approximates,  in  some  practically  useful  way,  a  map 
of  the  distribution  of  scattering  objects  in  the  x-y  plane. 

The  above  conditions  are  satisfied  if  x  and  y  are  taken  to 
be  time  delay  and  Doppler  frequency.  Determination  of  i]A(x,y) 
is  assumed  to  be  the  objective  of  the  radar. 

A  series  expansion  of  h(t;x,y)  is  obtained  in  terms  of  a 
set  of  functions,  $n(t),  which  are  orthonormal  over  the  re¬ 
ceiver  operating  time  interval  T,  and  a  set  of  functions 
0n(x,y),  which  are  orthonormal  over  a  finite  region  R  of 
the  x-y  plane.  If  the  only  a  priori  information  is  that 
i)r(x,y)  is  zero  outside  the  region  R,  then,  even  in  the  absence 
of  noise,  only  the  component  of  \|r(x,y)  which  is  representable 
as  a  linear  combination  of  the  0n(x,y)  can  be  deduced  from 
the  radar  return. 

A  smoothed  form  ^s(x,y)  of  \JKx,y)  is  defined  as  the  convolu¬ 
tion  of  \j)(x,y)  with  a  spike-like  smoothing  function.  In  the 
presence  of  additive,  stationary,  white,  Gaussian  noise,  a 
maximum  likelihood  estimate  of  (x ,y )  can  be  obtained  pro¬ 
vided  the  0n(x,y)  are  sufficiently  complete  to  permit  expand¬ 
ing  the  desired  smoothing  function  in  terms  of  the  0„(x,y). 

The  appropriate  receiver  is  not,  in  general,  a  matched  filter. 
The  cross  section  of  the  smoothing-function  spike  must  be 
chosen  to  compromise  between  loss  of  detail  in  i|rg(x,y)  and 
reduction  of  the  s ig na 1-t o-no i se  ratio. 

In  time  delay  and  Dopp ler-f requency  coordinates,  let  \JrT(x,y) 
be  the  component  of  ^(x,y)  that  contributes  to  the  radar 
return  during  the  receiver  operating  time  interval  T.  If 
it  is  known  a  priori  that  \jrT(x,2^k/T)  is  zero  for  x  outside 


i  i  i 


an  interval  D  and  for  2vk/T  outside  an  interval  B,  and  if 
the  transmitted  waveform  is  periodic  with  P  =  D,  then,  in 
the  absence  of  noise,  ^(x,2Trk/T)  can  be  deduced  exactly 
from  the  received  waveform  if,  and  only  if,  BD  <  2v  and 
all  harmonics  are  present  in  the  transmitted  waveform. 

If  the  smoothing  along  the  time-delay  axis  is  specified  by 
a  spike-like  function  a(x)  for  all  Doppler  frequencies,  and 
if  aCx)  is  a  realizable  autocorrelation  function,  then,  for 
a  given  transmitted  average  power,  the  output  s ignal-to-noise 
ratio  is  a  maximum  when  the  time  autocorrelation  function  of 
the  transmitted  waveform  is  a  periodic  string  of  pulses  of 
the  form  of  a(x).  The  receiver  weighting  function  appropriate 
for  such  a  waveform  resembles  closely  the  matched  filter. 
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CHAPTER  I 


INTRODUCTION 

A.  DELINEATION  OF  THE  AREA  OF  RESEARCH 

The  purpose  of  this  research  is  the  development  of  a  radar 
theory  applicable  to  dense  scatterer  distributions.  The 
scatterer  distribution  must  be  representable,  in  the  sense 
described  in  Chap.  II,  Sect.  A,  by  a  scatterer-density  function 
of  two  coordinates  x  and  y.  These  could  be  time  delay  and 
Doppler  frequency,  or  any  other  coordinate  pair  satisfying  the 
conditions  specified  in  Chap.  II.  Whether  or  not  other  such 
coordinate  pairs  exist  is  not  considered  in  this  report. 


B.  CRITIQUE  OF  EXISTING  MULTIPLE -TARGET  RADAR  THEORY 

A  fairly  complete  theory  of  detection  and  coordinate  estima- 

12  3  4  5  6* 

tion  is  available  for  the  case  of  a  single  point  scatterer.  ’ 

However,  in  most  radar  applications  of  interest,  there  is  more 
than  one  scatterer  present  at  any  one  time;  and  often  the  problem 
of  resolving  or  separating  any  one  scatterer  echo  from  other 
scatterer  echoes  with  arbitrary  coordinates  becomes  the  over¬ 
riding  consideration  in  radar  design.  Resolvability  (that  is, 
the  degree  to  which  scatterer  echoes  can  be  separated),  in  two 

* Superscripts  refer  to  numbered  items  in  the  Bibliography. 
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coordinates,  such  as  time  delay  and  Doppler  frequency,  is  an 
important  problem  in  multiple-scatterer  theory.  In  experiment¬ 
ing  with  different  transmitted  waveforms  ,  it  becomes  apparent 
that  good  resolvability  between  some  targets  may  be  achieved  at 
the  expense  of  resolvability  among  other  targets.  Often,  for  a 
given  transmitted  waveform,  the  resolvability  may  be  improved  at 
the  expense  of  detectability.  In  spite  of  these  observations, 
it  has,  so  far,  been  impossible  to  obtain  a  set  of  necessary  and 
sufficient  realizability  conditions  or  limiting  relationships 
which  prevent  the  simultaneous  achievement  of  arbitrary  resolva¬ 
bility  and  detectability  for  a  general  class  of  target  distribu¬ 
tions.  Furthermore,  a  systematic  procedure  for  the  synthesis  of 
the  receiver  and  the  transmitted  waveform  is  lacking. 

The  most  useful  guide  in  choosing  the  transmitted  waveform 

12  7 

has  been  Woodward's  uncertainty  function.  *  Siebert  has 
studied  the  properties  of  this  function  to  assist  the  synthesis 
of  the  transmitted  waveform  by  obtaining  a  test  of  the  realiza¬ 
bility  of  any  particular  form  of  the  uncertainty  function. 
Several  necessary  conditions  were  established,  but  the  lack  of 
a  useful  sufficient  condition  is  still  a  major  obstacle  in  the 
synthesis  of  the  transmitted  waveform. 


C.  THE  LIMITING  FORM  OF  THE  MULTIPLE -TARGET  PROBLEM 

It  appears  appropriate,  for  certain  radar  applications,  to 
consider  dense-scatterer  distributions  containing  scatterers  at 
all  points  throughout  a  region  of  the  plane  of  the  coordinates 
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x  and  y  in  which  the  distribution  is  to  be  mapped.  Under  these 
conditions  the  detection  of  any  one  scatterer  and  its  resolvabi¬ 
lity  from  all  the  other  scatterers  are  not  usually  of  direct 
interest  to  the  radar  user.  The  user  is  primarily  interested 
in  characteristic  patterns  formed  by  clusters  or  assemblies  of 
scatterers  and  their  over-all  distribution  in  the  x-y  plane. 

It  is  therefore  assumed  that  the  purpose  of  a  radar  system  in 
this  application  is  to  obtain  an  x-y  display  of  the  distribution 
of  scatterers.  Consequently,  the  appropriate  measure  of  perform¬ 
ance  is  the  fidelity  with  which  such  a  display  reproduces  the 
actual  scatterer  distribution. 

In  applying  existing  multiple-target  radar  theory  to  this 
dense-scatterer  problem,  the  following  difficulties  are 
encountered : 

1.  The  theory  is  incomplete;  it  lacks  limiting 
relationships,  realizability  conditions,  and 
a  systematic  synthesis  procedure. 

2.  The  measures  of  performance,  that  is  detect¬ 
ability  and  resolvability,  are  not  directly 
related  to  the  interests  of  the  radar  user. 

As  a  result,  a  new  approach  to  the  dense-scatterer  problem 
has  been  undertaken. 


D.  THE  SCOPE  OF  THE  PROPOSED  THEORY 

An  appropriate  scatterer-density  model  representing  the 
physical  scatterer  distribution  is  developed.  The  determination 
of  this  model  can  then  be  regarded  as  the  objective  of  the  radar 
system  in  dense-scatterer  applications.  Limiting  relations 
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showing  the  extent  to  which  the  scatterer-density  model  may  be 
deduced  from  the  radar  return  are  obtained  for  the  noise-free 
and  the  noisy  cases.  A  systematic  synthesis  procedure  for  the 
appropriate  receiver  is  developed.  Also,  the  synthesis  of  trans¬ 
mitted  waveforms  is  carried  out  for  the  case  in  which  x  is  time 
delay  and  y  is  Doppler  frequency. 

The  proposed  theory  is  portinent  to  applications  involving 

e 

dense-scatterer  distributions.  The  single-scatterer  case  has 
been  thoroughly  covered  by  existing  radar  theory.  Between  the 
single-scatterer  and  the  dense-scatterer  distributions,  there 
lies  a  large  intermediate  class  of  distributions  which  may  be 
approximated  by  a  finite  number  of  discrete  point  scatterers  , 
whose  coordinates  are  not  known  £i  priori .  With  respect  to  these 
distributions,  the  dense-scatterer  theory,  just  as  the  single- 
scatterer  theory,  may  prove  useful  as  a  limiting  form  of  the 
actual  problem. 


CHAPTER  II 


GENERAL  THEORY 


A.  THE  RADAR  MODEL 

The  bandwidths  of  all  pertinent  radar  waveforms  (the 
transmitted  waveform,  the  return  from  a  point  scatterer,  and  the 
return  from  an  arbitrary  distribution  of  scatterers)  are  assumed 
to  be  less  than  twice  the  carrier  frequency  so  that  each  waveform 
may  be  represented  as  a  complex  video  signal  with  a  one-one 
relation  obtaining  between  waveforms  and  their  representations. 
Such  a  representation  does  not  contain  the  carrier;  only  the 
amplitude  and  phase  modulation  are  preserved.  Suitable  linear 
combinations  of  the  real  and  imaginary  parts  of  the  complex 
video  yield  in-phase  and  quadrature  components  of  actual  video 
signals  encountered  in  coherent  radar. 

The  complex  video  waveform  representing  the  return  from  a 
point  scatterer  is  designated  by 

Return  from  a  point  scatterer  =<^h(t;x,y)  (2-1) 

where  h(t;x,y)  depends  on  the  radar  characteristics 
(transmitted  waveform,  antenna  pattern,  and  antenna  scan),  t  is 
time,  x  and  y  are  the  coordinates  of  the  scatterer,  and  ¥  is  a 
complex  amplitude  representing  the  magnitude  and  phase  of  the 
return. 
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The  return  from  a  number  of  point  scatterers  so  located 
that  the  return  from  each  is  independent  of  the  presence  of  the 
others  is  given  by  a  linear  superposition  of  returns  of  the 
form  of  Bq.  2-1;  namely, 


Return  from  a  set  of 
non-interferring  point  scatterers 


where^^,  xk  ,  and  y^  are  the  complex  amplitude  and  the  coordi* 
nates  of  k*''1  scatterer.  For  dense-scatterer  distributions 
containing  scatterers  at  all  points  of  the  x-y  plane,  the  analog 
of  Eq.  2-2  is: 


00 

s(t)  =  J  J  NMx.y)  h(t  ;  x  ,y)  dx  dy  (2-3) 

-  00 

where  s(t)  is  the  radar  return  andY  (x,y )  is  the  density  per 
unit  area  of  the  x-y  plane,  of  the  sources  of  backscatter  at 
a  point  (x,y).  Discrete  point  scatterers  may  be  accounted  for 
by  two-dimensional  impulses  in Y(x,y)  at  the  appropriate  points 
in  the  x-y  plane.  Consequently,  Eqs.  2-1  and  2-2  may  be  con¬ 
sidered  special  cases  of  Eq .  2-3.  Thus  for  noninteracting 
scatterers,  a  display  of  jY(x,y)|  is  a  map  of  the  distribution 
of  physical  objects  in  x-y  coordinates  — nonzero  and  zero  values 
of  |S/’(x,y)|  indicate,  respectively,  the  presence  or  absence  of 
objects  at  the  corresponding  regions  of  the  x-y  plane. 

Configurations  of  objects  of  practical  interest  often 
consist  of  interacting  scatterers  which,  among  other  things,  may 
give  rise  to  destructive  interference  and  fading  of  the  target, 
the  appearance  of  false  targets  due  to  multiple  reflection,  and 
loss  of  targets  due  to  shadowing  by  objects  at  closer  range. 
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Despite  this,  it  may  still  be  possible  that  in  an  appropriate 
set  of  coordinates,  Eq.  2-3  is  valid  for  a  function  h(t;x,y) 
determined  by  the  radar  and  a  function  SMx ,y)  that  depends  on 
the  scattering  objects  and  that,  under  the  proper  conditions 
(such  as  those  discussed  in  Seot .  A  of  Chap.  Ill),  can  be  made 
independent  of  the  radar.  It  is  shown  in  Sect.  A  of  Chap, 

III  that  time  delay  and  Doppler  frequency  are  an  appropriate  set 
of  coordinates  which  meet  the  above  requirements. 

In  consequence  ,  we  consider  henceforth  some  coordinate  pair 
x,y  for  which  we  are  able  to  define  a  scatterer-density  function 
haring  the  following  important  properties: 

1.  NMx.y)  is  independent  of  the  radar 

2.  the  received  waveform  for  a  radar  characterised 
by  any  function  h(t;x,y)  is  given  by  Eq.  2-3 

3.  Y(x ,y)  characterises  a  class  of  actual  scatterer 

configurations  having  members  which  are  indis¬ 
tinguishable  by  a  radar  characterised  by  any 
h(t;x,y)  and  using  any  receiver  processing 

4.  |^(xfy)|  may  sometimes  be  interpreted,  at  least 
approximately,  as  a  distribution  of  objects 
throughout  the  x-y  plane 

In  view  of  the  above  properties  and  because,  given  S' ( x  ,  y )  ,  |Y(x,y) 
is  available  for  whatever  interpretation  as  a  map  of  che  distribu¬ 
tion  of  objects  in  the  x-y  plane  may  be  possible,  the  desired 
radar  output  is  assumed  in  this  report  to  be  an  x-y  display  that 
reproduces  NMx  ,y)  within  some  prescribed  fidelity  criterion. 

We  have,  in  effect,  subdivided  the  dense-scatterer  problem 


into  three  parts: 
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1.  the  determination  of  a  suitable  pair  of  coor¬ 
dinates  and  a  corresponding  f unction  Y(x  ,y ) 

2.  the  interpretation  of|Y(x,y)|  and  its  rela¬ 
tionship  to  the  distribution  of  objects 

3.  the  establishment  of  conditions  under  which 
the  desired  radar  output  is  achievable  and 
the  determination  of  the  appropriate  trans¬ 
mitted  waveform,  antenna  pattern  and  scan, 
and  receiver  processing  functions  for  obtain¬ 
ing  the  desired  output. 

Henceforth,  aside  from  Sect.  A  in  Chap.  Ill,  our  only  concern 
will  be  part  3  of  the  dense-scatterer  problem. 


B.  PERFORMANCE  LIMITATIONS  IN  THE  ABSENCE  OF  NOISE 

In  the  following  we  shall  determine  the  extent  to  which  the 
scatterer  density  can  be  deduced  from  the  knowledge  of  the  noise 
free  received  and  transmitted  waveforms  if  the  only  available 
JL  Priori  information  is  that  Y(x,y)  vanishes  outside  some  finite 
region  R  of  the  x-y  plane.  In  other  words,  we  shall  inquire 
into  the  possibility  of  solving  for  Y(x ,y)  from  Eq .  2-3  when 
s(t)  and  the  transmitted'  waveform,  and  consequently  h(t;x,y), 
are  known  and  the  region  of  integration  is  restricted  to  a 
finite  region  R. 

First,  two  basic  sets  of  orthonorraal  functions,  one  in  t 
and  one  in  x  and  y,  will  be  obtained.  The  radar  return  s(t), 
the  return  from  a  point  scatterer  h(t;x,y),  and  the  scatterer 
density Y ( x ,y)  will  then  be  represented  in  terms  of  the  two 

basio  sets.  Vlth  the  help  of  these  representations,  the 
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solution  to  the  problem  posed  above  will  become  obvious. 
Furthermore,  the  two  basic  sets  of  orthonormal  functions  will 
prove  very  useful  in  obtaining  ultimate  performance  limitations 
in  the  presenoe  of  noise  and  in  developing  a  synthesis  procedure. 
Ve  begin  by  defining  a  function  K(t,t')  as  follows: 


K(t  ,t  •  )  ■  JJ  h(t ;x ,y )  h*(t';x,y)  dx  dy  (2-4) 

R 

where  the  asterisk  denotes  the  complex  conjugate.  It  is  assumed 
that  |h(t;x,y)|  is  bounded  and  that  the  region  R  is  finite.  Ve 
note  that 

K(t,t'  )  =  K*  (t 1  ,t) 

and  that,  as  a  oonsequence  of  the  above  assumptions. 


/  /  lK<‘  ,  t ' ) 1  dt  dt' 

T  T 

is  finite  for  any  finite  T.  Thus,  K(t,t')  is  an  (integrable 

square),  Hermitian  kernel.  ,7  The  corresponding  characteristic 

functions  cp  (t)  and  characteristic  values  A  satisfy  the  follov- 
n  n  * 

ing  equations:* 

f  K(t,t«)  ?n(t')  dt»  -  An  cpn(t) 

T 

f  < Pm(t)  ,Pn*(*) 

T 

A 

n 


■ 

\  *„<*> 

n-1,2,3  , . . . ; 

t  inside 

(2-5) 

fi 

if 

m=n 

m 

{o 

if 

m^n 

(2-6) 

0 

n*l ,2  ,3. . . 

♦Most  treatments  of  integral  equations  consider  only  real 
symmetric  kernels  except  that  the  statement  is  made  that  the 
results  can  be  readily  extended  to  Hermitian  kernels  as  is  done 
here. 
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The  interval  of  integration  T  is  the  observation  interval  which 

is  a  fixed  period  of  time  for  whioh  the  received  waveform  s(t) 

(see  Eq.  2-3)  is  available  at  the  radar  receiver. 

Knowing  the  characteristic  functions  ?m(t),  we  generate  a 

set  of  funetions  6  (x,y)  as  follows: 

m 

J  h*(t';x,y)  dt*  ^  pm  eo(x,y)  (2-7) 

T 

where  the  M>m  are  normalising  constants  such  that 

*  .  2  • 

JJ  I6m*x,yH  dx  dy  -  1  (2-8) 

R 

Since  the  normalization  can  always  be  accomplished  with  a  posi¬ 
tive  real  constant,  we  shall  take  the  nn  to  be  positive  and  real. 
As  a  result  of  Eq.  2-7,  we  have 


J  h(t  ;x  ,y )  cpn*(t)  dt  “  9n*(x,y)  (2-9) 

T 

Multiplying  Eq.  2-7  by  Eq.  2-9,  integrating  with  respect  to  x 
and  y,  and  interchanging  the  order  of  integration  on  the  left- 
hand  side,  we  obtain 


J  H  ;x'y)h*  (tf;*.y)dxdy  <Pn*(t  )dtdt 1  = 

T  T  L  r  J 

*V*n  //  0m(x'y)  en*{x’y)dxdy 
R  (2-10); 

Substituting  Eq .  2-4  for  the  bracketed  term  on  the  left  -hand 
side  and  making  use  of  Eqs.  2-5,  2-6,  and  2-8,  we  find  that 


//  emU.y)  VU.y)  dx  dy  "  {  0 


if  m  =  n 


if  m  f  n 


(2-11) 
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and  that 


p  2  =  X 

*n  n 


(2-12) 


Equation  2-12  indicates  that  K(t,t')  is  a  non-negative  kernel 
(recall  that  the  ^n's  are,  positive  and  real),  which  fact  oan 
also  ba  ascertained  by  noting  that 

//  /h(t;x,y)  u(t)  dt I  dx  dy  *  0 
R  T 

which  with  the  help  of  Eq.  2-4  becomes 


J  /u(t>  K(t ,t 1 )  u*  (t '  )  dt  dt 1  *  0 


(2-12a ) 


for  any  u(t ) . 

It  is  useful  to  note  that  the  functions  cpn(t)  can  be  gene¬ 
rated  from  the  functions  6n(x,y)  by  first  multiplying  both  sides 
of  Eq.  2-7  by  h(t;x,y),  integrating  with  respect  to  x  and  y,  and 
interchanging  the  order  of  integration,  so  that 

M-n  J^htt  ;x  ,y  )0n  (x  ,y  )dxdy  *  J  JJ  h(t ; x ,y )h* (t ' ;x ,y )dxdy  cpn(t')dt' 
R  T  LR 

Next,  by  substituting  K(t,t')  as  given  by  Eq.  2-4  for  the  brack¬ 
eted  term  on  the  right-hand  side  and  by  using  Eqs.  2-5  and  2-12, 
we  obtain 


J J'hlt  ;x,y)  0n(x,  y)  dx  dy  *  un  9n< 


(2-13) 


In  order  to  obtain  a  representation  for  h(t;x,y)  in  terms 

of  the  ®  (t)  and  6  (x,y),  oonsider  the  following  expression 
n  in 
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en  *  fff\h(tix,y)  -  ^  ^(t)  ®k*(x,y)|  dxdydt  (2-l4) 

T  ft  k=l 

With  the  help  of  Eqs.  2-4,  2-13,  and  2-12,  wo  find  that 

n 

«n  =  /K(t  ,t )  dt  -  ^  \  (2-15  ) 

T  k=l 

It  will  be  assumed  that  the  radar  return  from  a  point  scatterer , 
that  is,  h(t;x,y),  is  a  continuous  function  of  time  uniformly  in 
x  and  y  over  R  —  an  assumption  which,  from  a  physical  standpoint, 
does  not  restrict  the  generality  of  the  analysis.  Consequently, 
K(t,t')  as  defined  by  Eq.  2-4  must  also  be  continuous.  Further¬ 
more,  since  K(t,t')  is  non-negative  (see  Eq.  2-12a),  we  can  use 
a  corollary  of  Mercer's  theorem  (see  reference  9»  p.  127),  which 
states 


oo 

)  Ak  *  fxit.t)  dt 


(2-16) 


k=l  T 

Therefore,  as  n  approaches  infinity,  Eq,  2-15  becomes 

4  im  e  =0 
-  n 
n  — ►  oo 

which  when  applied  to  Eq.  2-l4  yields 

n 

l  im  /// 1  »(*!*./>  -  YL  ^k  *Pk ^  |d*  dy  dt  »  0 

n  *■  °°  T  R  k=l 

Thus,  we  obtain  the  representation 

n 


h( 


tjx.y)  =  4 . i.m.  y  ^  <Pk  ( t )  6k*(x,y)  (Z~17) 


which  is  valid  regardless  of  the  completeness  or  incompleteness 
of  the  sets  9k(t)  and  0k(x,y). 
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Next,  we  seek  the  appropriate  series  representations  for 
the  radar  return  s(t).  Multiplying  both  sides  of  Eq .  2-3  by 
cpk*(t)  integrating  over  T,  interchanging  the  order  of  integra¬ 
tion  on  the  right-hand  side,  and  using  Eq.  2-9,  we  find  that 

•y.  ■  ek*U,y)  dx  dy  (2-18) 

R 

where 

sk  =  J s(t)  9k*(t)  dt  (2-19) 

T 

Multiplying  both  sides  of  Eq.  2-18  by  cpk  ( t )  and  summing  from 
k*l  to  k=n ,  we  obtain 
n  n 

y  .  ak  *  y  Hk  yk(t)  ek*(x.y)dxdy  (2-20) 

k=l  k=l  R 

Subtracting  Eq.  2-20  from  Eq.  2-3  and  re-arranging  terms  gives 
n  n 

s(t)  ~y  ■  sk^k  ^  ^  =  jj[^(t  ;x,y)  -/  M>k<Pk(t  )6k*(x,y)]  (|/(x,y)dxdy 
k= 1  R  k=l 

By  applying  the  Schwars  inequality,  we  obtain 

i«(t)  .y  ^iti|  * 

k=l 


2 

/T|h(t;x,y)  -  ^«Pk(t)ek*(x,y)|2  dxdy  ff\^x>y)\  d*dy 

R  k= 1  R  (2-21) 

With  the  help  of  Eqs.  2-4,  2-12,  and  2-13,  the  first  faotor  on 

the  right-hand  side  may  be  rewritten  so  that 

n  2  n  2 

(t) -Y2  sk9k<t)|  *[K(t,t)  -r.  Vk(t)cpk*(t)]  Jf  |^(  x,y)J  dxdy 

R  (2-21a) 


k*l 


k«l 
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Since  K(t,t')  is  continuous  and  non-negative  (see  discussion 
below  Eq.  2-15)  »  we  may,  in  proceeding  to  the  limit  as  n 
approaches  infinity,  apply  Mercer's  theorem  to  the  bracketed 
expression  in  Eq.  2-21a.  Thus,  assuming  that  \f/(x,y)  is  an 
integrable-square  function*  over  the  region  R,  we  obtain 


so  that 


oo 

Bit)  •  sk  Cpk(t)  (2-21b) 

k=l 
9 

converges  almost  uniformly.  Note  that  Eq .  2-21b  is  valid 
regardless  of  the  completeness  or  incompleteness  of  the  cpk(t). 

We  are  now  in  a  position  to  inquire  into  the  possibility 
of  solving  for  \p(x,y)  from  Eq.  2-3.  For  this  purpose  we  define 
the  coefficients  ^  as  follows: 


®  Jf  ^(x.y)  0k*(*.y)  d*  dy  (2-22) 

R 

where  we  have  used  Eq.  2-18.  From  the  theory  of  orthonormal 
expansions^,  it  is  known  that  \j/(x, y)  may  be  decomposed  into  two 
orthogonal  components  as  follows: 

t(*.y)  ■  fy|(**y)  +  <^(x,y)  (2-23) 

-such  that 


♦  Such  a  restriction  excludes  point  targets  which  produce  impulses 
in  \f/(x,y).  It  does  not  exclude  approximations  to  point  targets 
that  produce  large -magnitude  narrow-base  spikes  in  ^(x,y). 
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ff  |{(*,y)  fj*(x>y)  dx  dy  ■  0 
R 


ekU.y) 


//|4'l|(*.y)|  d*  “x  *  XI  Kl 

R  k-1 


(2-23*  ) 


(2-23b) 


(2 -23c ) 


Jff(x,y)  6k*(x,y)  dxdy  -  Jf^\\^x,y^  ®k**x,y*  dxdy  * 

R  R  (2-23d) 

//*±<x  #y)  S^^tx.y)  dxdy  *  0  for  each  k  (2-23e) 

R 

and 

II  +  11  *  II  t||  INI  till  (2-23f) 

where  the  double  vertical  bars  are  used  to  denote  the  norm  or 
integrated  square  of  the  function  shown  between  the  bars.  If 
the  system  of  functions  6k(x,y)  is  complete  over  the  region.  R, 
then  \f/ jjx.y)  ■  0. 

By  substituting  Eq.  2-23  into  Eq.  2-18  and  making  use  of 
Eq.  2-23e,  we  find  that  \ f/jjx,y)  makes  no  contribution  to  any  of 
the  sk*  Thus,  from  Eq.  2 -21b,  it  follows  that 

JJ h(t;x,y)  ifrjJ x  , y )  dx  dy  “  0  (2-24) 

R 

Consequently,  the  component  ^(x,y),  which  will  be  referred  to 
as  the  ambiguous  component  of  the  scatterer  density,  does  not 
contribute  anything  to  the  radar  return;  and  the  most  that  we 
can  hope  to  determine  from  s(t)  is  the  unambiguous  component 
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(x,y).  The  latter  can  indeed  be  determined  by  a  radar 
receiver  which  operates  on  s(t)  with  a  weighting  funotion 
Vn(x,y;t)  defined  by 

n 

V*,y;t)  ®k(x,y)  (2-25) 

k^T  k 

so  that  the  result  of  the  operation  is 

n 

/wnU.y;t)  s(t)  dt  =  y  ek(x,y)  (2-26) 

T  k  =  l 

where  we  have  used  Eq.  2-3,  2-9,  and  2-22.  In  the  limit  as  n 
approaches  infinity,  we  obtain  the  representation  for  >K||(x,y) 
given  by  Eq.  2 -23b. 

Recall  that  K(t,t')  is  a  non-negative  L 2  kernel  so  that  its 

characteristic  values  *n ,  and  hence  the  un  (see  Eq.  2-12),  may 

8  9 

be  arranged  in  a  sequence  of  nonincreasing  values  with 
increasing  n.  Consequently,  if  there  is  an  infinite  number  of 
Un's,  the  weighting  function  given  by  Eq.  2-25  does  not  tend  to 
a  limit  as  n — »  oo .  However,  the  result  of  operating  with  this 
weighting  function  (that  is,  Eq.  2-26)  does  have  a  limit  in  the 
mean  by  Eq .  2-23b. 

In  the  foregoing  we  have  obtained  a  representation  of 

h(t;x,y)  in  terms  of  its  characteristic  functions  cp  (t)  and 

n 

8n(x,y)  and  its  characteristic  values  nn.  The  scatterer  density 
has  been  decomposed  into  two  orthogonal  components  Y||(x,y)  and 

HMx.y)  such  that  the  first  is  representable  as  a  linear  combina- 

•  ~ 

tion  of  the  characteristic  functions  of  h(t;x,y);  whereas,  the 
second  is  not.  Thus,  the  radar  return,  which  was  shown  to  con¬ 
tain  contributions  only  from  t-he  characteristic  functions  of 
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h(t;x,y)  carries  no  information  about  YjJx ,y) .  As  a  result, 
only  Y( | ( x , y )  can  be  deduced  from  the  radar  return.  For  this 
reason  Y||(x,y)  and  YjJx.y)  are  referred  to,  respectively,  as 
the  unambiguous  and  the  ambiguous  component's  of  the  scatterer 
density.  If  h(t;x,y)  and  the  region  R  are  such  that  the 
characteristic  functions  0n(x,y‘)  are  complete,  then,  in  the 
absence  of  noise,  a  series  representation  for  the  total 
scatterer  density Y (x ,y)  can  be  deduced  from  the  radar  return. 

The  completeness  condition  may  be  restated  in  a  number  of 
8  9 

equivalent  forms.  *  However,  there  is  no  direct  mathematical 
test  for  explicitly  establishing  the  completeness  or  incomplete¬ 
ness  of  the  8  (x,y)  for  any  h(t;x,y).  Thus,  the  problem  of 
deducing  Y(x ,y )  from  s(t)  has  been  reduced  to  a  well-known 
mathematical  problem. 

The  indexing  of  the  characteristic  functions  0k(x,y)  is 
arranged  in  order  of  nonincreasing  pk  or  ^k  which,  in  general, 
can  be  assumed  to  correspond  to  increasing  "wiggliness"  of  the 
8k(x,y)  as  a  function  of  x  and/or  y.  Thus,  high-order  (high-k) 
characteristic  functions  represent  the  sharp  features  of  Y||(x,y) 
whereas,  low-order  (low-k)  characteristic  functions  represent 
the  d-c  -  like,  or  smooth,  features  of  Y||(x,y). 

The  incompleteness  of  a  set  of  characteristic  functions 
0k(x,y)  may  be  regarded  as  a  deficiency  of  orthogonal  functions 
required  for  the  representation  of  Yj^(x,y).  The  set  of  miss-ing 
functions  can  also  be  classified  into  high-  and  low-order 
functions.  Deficiencies  of  the  first  class  result  in  an  error 
in  the  sharpness  of  detail  on  the  radar  display  (which  can  only 
show  Y||(x  .y)  ).  and  deficiences  of  the  second  class  result  in 
an  error  in  the  smooth  features  on  the  display. 
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It  is  clear  from  Eq.  2-23  that  ignorance  of  (//jjx, y)  will 
also  make  it  impossible  to  obtain  a  display  of  jl//(x,y)J.  Quali¬ 
tatively,  the  deficiencies  in  a  display  of  |^||(x,y)|  are  the 
same  as  the  deficiencies  in  ^||(x.y)  discussed  above. 

C.  SMOOTHING  AND  ESTIMATION  OF  THE  SCATTERER  DENSITY  IN  THE 
'  PRESENCE  OF  NOISE 

In  the  presence  of  noise,  which  is  assumed  to  be  Gaussian, 
stationary,  white,  and  additive,  it  becomes  impossible  to  make 
an  exact  determination  of  ^||(x,y)  or  |^||(x,y)|  ;  .  much  less, 

of  \J/(x,y)  or  |i//(x,y)|.  It  will  be  shown  that,  under  these 
conditions,  only  ah  estimate  of  a  smoothed  form  of  <//j|(x,y)  can 
be  obtained.  Such  behavior  is,  to  a  certain  extent,  analogous 
to  that  encountered  in  the  problem  of  extracting  an  unknown 
signal  s(t)  from  a  background  of  white,  Gaussian,  additive  noise. 
In  order  to  obtain  an  undistorted  replica  of  s(t),  an  infinite 
bandwidth  is  required  which,  in  turn,  results  in  infinite  noise 
power.  Therefore,  a  certain  degree  of  smoothing  and  the  result¬ 
ing  signal  distortion  is  accepted  in  order  to  reduce  the  noise 
power.  The  optimum  smoothing  operation  (that  is,  one  that  mini¬ 
mizes  some  measure  of  the  total  error  caused  by  signal  distortion 

/ 

and  noise)  depends  on  the  characteristics  of  s(t)  and  can  be 
determined  only  from  ja  priori  knowledge  about  s(t).  Since  such 
information  is  not  available  for  ^(x,y),  we  shall  consider  the 
effect  of  a  class  of  smoothing  operations  upon  both  signal,  in 
our  oase  ^(x,y),  and  noise. 
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One  example  of  a  linear  smoothing  operation  is  given  by: 

Ys(x,y)  -  Jf  S(x-x*  ,y-y' )  Y(x' ,y' )  dx'  dy«  (2-2?) 

R 

where  Vg (x ,y)  is  the  result  of  smoothing  and  S(x,y)  has  the  form 
of  the  right  cylinder  illustrated  in  Fig.  la.  The  cross-section 
area  of  the  cylinder  is  O,  the  height  is  1  jo ,  and  the  volume  is 
one.  For  such  an  S(x ,y) ,  Y  (x ,y )  in  Eq.  2-2?  is  the  average  of 
Y(x  ,.y)  over  the  area  cr  . 

More  generally  |s(x,y)|  may  have  any  spike-like  form,  such 
as  the  function  illustrated  in  Fig.  lb.  The  result  of  smoothing, 
in  this  more  general  case,  can  be  thought  of  as  a  weighted 
average  ofY(x,y),  if  the  volume  under  S(x,y)  is  set  equal  to  one 

oo 

//  S(x,y)  dx  dy  *  1  (2-28) 

-  oo 

It  is  convenient  to  define  the  effective  averaging  area  of  a 
general  spike-like  S(x,y)  as  the  cross  section  O'  of  a  cylinder 
that  is  equivalent  to  the  spike  in  the  sense  that  the  volumes 
under  the  squared  magnitude  of  the  two  functions  are  equal.  Thus 

cr  -  - i -  (2-29) 

00  2 
j J  | S(x ,y ) |  dx  dy 

-  oo 

Of  particular  importance  in  the  present  context  is  the  class  of 
smoothing  functions  for  which  (a)  S(x,y)  falls  off  rapidly ‘for 
x  and  y  .outside  an  area  o  about  x  *  0  ,  y=  0;  and  (b)  O  is  much 
smaller  than  the  area  of  the  region  R.  For  such  a  smoothing 
function,  the  radar  dleplay .  oontains  significant  detail.. 
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If  we  regard  S(x-x' ,y-y* )  as  a  function  of  X*  and  y'  with 
jr  and  y  as  parameters we  oan  write 

S(x-x'.y-y')  ■  y~Pn(x.y)  «„*(*'  ,y' )  +  Sjx  ,y  ;x*  ,y ' )  (2-30) 

n 

where 

Pn(x,y)  =  JJ*  S( x-x'  ,y-y'  )  Sn  ( x '  ,y'  )  dx'dy'  (2-30a) 

R 

and  Sj^(x,y;x'  ,y‘ ).  is  the  component  of  $(x-x' ,y-y' )  orthogonal  to 
all  the  functions  8  (x',y'),. 

Jf Sj(x,y;x' ,y»)  8n(x',y')  dx'dy'  -  0  (2-30b) 

R 

Substitution  of  Eq.  2-30  into  Eq.  2-27  and  use  of  Eqs.  2-22,  2-23, 
2-23b,  and  2~30b  yield 

N'Jx.y)  “^YkPk(x,y)  +  JJsjU.yjx*  ,y«)  Yj_(x'  ,y«)  dx'dy' 

k  R  (2-31) 

Since  only  theYk's  can  be  deduced  from  the  received  waveform 
(even  under  noise-free  conditions),  the  last  term  in  Eq.  2-31 
constitutes  an  ambiguous  component  of  Ys(x,y).  One  means  of 
obtaining  unambiguous  determination  of  Y  (x,y)'  in  the  noise-free 
case,  .  therefore,  is  to  select  a  smoothing  function  and  trans¬ 
mitted  waveform  in  such  relationship  to  one  another  that 
S^(x  ,y ;x' ,y' )  is  sero.  In  these  circumstances, 

Ys(x.y)  =  ^Yk  pk(*.y)  (2-3ia) 

k 

The  radar  return  from  which  Y  (x,y)  is  to  be  estimated  is 

S 

r(t)  *  s(t)  +  n(t)  (2-32) 


given  by: 
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where  s{t)  is  the  signal  component,  given  by  Eq.  2-3,  and  n(t) 
is  the  noise  component  in  tho  complex  video  representation.^ 
Multiplying  both  sides  of  Eq.  2-32  by  cp^ * ( "t )  and  integrating 
with  respect  to  t  over  the  interval  T  gives 

rk  "  *k  +  nk  k  «  l,  2,  3,  ...  (2-33) 

where 

rk  *  Jr  ( t )  cpk  *  ( t )  dt  (2  -33a  ) 

T 

s^  is  given  by  Eq.  2-19,  and 

nk  =  dt  (2-33b) 

T 

With  the  help  of  Eq.  2-22,  Eq .  2-33  may  be  rewritten  as  follows: 

rk  "  “kYk  +  nk 

Equation  2-33b  is  to  be  interpreted  as  a  stochastic  integral10; 
and  nk y  a  random  variable. 

Since  the  s^1 s  completely  specify  the  signal  component  s(t) 

(see  Eq .  2-19),  the  r^ ' s  may  be  used  as  the  "observables"  or 

3  4 

"observable  coordinates."  *  Thus,  the  problem  of  estimating 
Ys  ( x ,  y )  from  r(t)  may  be  restated  as  a  problem  of  estimating 
N^(x,y)  from  the  r^'  s. 

It  was  assumed  that  the  noise  n(t)  is  stationary,  white, 
Gaussian  noise  of  zero  mean;  let  its  spectral  density  be  Nq . 

With  the  help  of  Eq.  2-33b,  we  can  evaluate  the  following 


averages : 
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//nltT  'WJ 

T  T 


<Pk(t)  dt  dt1 


■  //» 0 

T  T 


6  C t •  )  9k(t)  <pA*(t»)  dt  dt' 


nk  n£ *  =  N0  /  ‘Pic***  <PA * < *  >  dt 


(2-35) 


where  the  horizontal  bar  denotes  an  ensemble  average.  Thus, 
using  Eq.  2-6,  we  obtain 

if  k  =  l 
if  k  it 


n*n‘*  =  {o° 


(2-36) 


Since  linear  operations  on  Gaussian  random  processes  result 


in  Gaussian  random  processes  or  variables, 


4,10 


all  the  nk's  are 


jointly  Gaussian  random  variables  with  zero  mean.  Thus,  for  a 
given  set  of  ^'s,  the  r^'  s  as  given  by  Eq.  2-34  are  jointly 
Gaussian  random  variables  whose  joint  probability  density  will 
be  denoted  by 

Probability  density  -  p(ri  »r2  »  •  *  ./^  ,^2  .^3  •  *  •  • )  (2”37) 


From  Eqs .  2-34  and  2-37  and  from  the  fact  that  the  average  value 
of  n^  is  zero,  we  conclude  that  the  maximum-likelihood  estimate 
^ce  of  ^k  ls  Siven  by 


(2-38) 


If  S(x-x'  ,y-y'  )  is  nonsingular,  that  is*  if,  given  l^g(x,y)  , 
Eq.  2-30  is  satisfied  by  a  unique  function  v|/(x, y);  and 
furthermore,  if  Sj^(x,y;x'  , y 1  )  is  zero,  then  a  one-to-one 
correspondence  holds  between  sets  of  and  functions 
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Ys(x.y) ,  as  given  by  Eq.  2-31a.  Under  these  circumstances  the 

maximum-1 ikelihood  estimate  of  Ys(x,y)  is  formally  given  by 

oo 

■  X]Yke  Pk<x-y>  <2-3S>) 

k=l 


where  Yge  (x  ,y)  is  the  estimate  .of  Yg(x,y).  Using  Eqs.  2-38,  2-34, 
and  2— 31a  in  Eq,  2-39  and  redefining  terms,  we  obtain 

^e^.y)  =  Ys(x,y)  +  n (x ,y )  (2-40) 

where  n(x,y)  is  given  by 


n (x ,y ) 


E 

\^  =  i 


(x  ,y) 


( 2 -40a ) 


Thus,  according  to  Eq.  2-40,  the  final  output  on  the  radar  dis¬ 
play  consists  of  two  components,  one  of  which  is  the  smoothed 
scatterer  density  and  the  other  is  the  noise.  From  Eqs.  2-40, 
2-40a,  2 -30a  ,  and  2-27,  it  follows  that 


Yse(x.y)  =  //s(x-x\y-y’)YU'  ,y*  )  dx'dy* 
R 


+  E 


Jfsix-x 

-  R 


.y-y' ) 


0k(x' ,y 1 )  .  dx  'dy' 


-1(2-41) 


If  S  ( x-x 1  ,  y-y11  )  is  s  ingular--that  is,  if  knowledge  of  Y  (x,y) 

s 

doesnot  det  ermine  Y(  x  ,y)  uniquely,  but  still  Sj_(  x  ,  y  ;  x '  ,y 1  )  is 
zero-r-the  set  ofYk's  cannot  be  deduced  from  Ys  ( x  ,  y )  .  Observe, 
however,  that:  (l)  an  estimate  Y  (x  ,y)  is  determined  by  Eq .  2-4l, 
whether  or  not  S(x-x' ,y-y' )  is  singular;  and  (2)  as  the  input 
noise  decreases  to  zero,  the  n.  's  tend  to  zero  and  Y  (x,y) 
approaches Ys ( x.y ) ,  whether  S(x-x'  ,y-y'  )  is  singular  or  not.  We 


( 

I 

i 

4 

i 

i 

l 

\ 

-f 
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therefore  consider  Eq.  2-41  to  yield  a  reasonable  estimate  of 
^8(x,y )  in  the  singular  case. 

If  Sj_(  x  ,y ;  x '  ,y  ’  )  is  not  zero,  the  set  of  0^  ( x 1  ,  y '  )  is  not 
complete;  and  for  some  »j/{x'  ,y'  ),  yf/jlx'  ,y'  )  not  only  is  not  aero 
but  also  is  not  orthogonal  to  S^(x,y;x' ,y' )  for  some  (x,y).  By 
Eq.  2-31,  therefore,  an  ambiguous  component  of.^g(x,y)  exists 
even  in  the  noise-free  case.  We ,  therefore,  conclude  that  a 
meaningful  estimate  {given  by  Eq.  2-41)  can  be  obtained  if, and 
only  if,  the  set  of  0^(x,y),s  is  sufficiently  complete  so  that 
S(x-x' ,y-y' ) ,  regardless  of  whether  it  is  singular  or  not,  is 
representable  in  the  form  of  Eq.  2-30b. 

The  mean-square  value  of  the  noise  at  any  point  on  the 
display  is  given  by 


_ _______ __  00 

|n(x,y)|2  =  Nq 


Pk(x,y) 


(2-42) 


k=l  '“k 

where  we  have  made  use  of  Eqs.  2-30a  and  2-36  in  Eq.  2-40a.  If 
we  average  the  mean-square  noise  over  the  area  of  the  display, 
that  is,  the  region  R,  we  obtain  the  average  mean-square  noise 


___________  00 

nr  "  r  fj  l« *y > 1 2  *  NQ  ^ 


(2-43) 


k=l 


he’ 


2  2 
where  pk  is  the  average  value  of  |pk(x,y)|  ;  namely, 

pk2  “  R  //lpk(x‘y)|2  dxdy  (2-43a) 

R 

In  applications  in  which  |vj/(x,y)|  is  a  meaningful  representa 
tion  of  the  scatterer  distribution  (see  Sect.  A  of  this  chapter), 
the  radar  user  may  be  primarily  interested  in  the  magnitude  of 


2  6 


Y  (x,y).  Unfortunately  the  relation  between  Y_(x.y) 


and  the 


set  of Yk' s  is  not  one-to-one  so  that  a  maximum-likelihood  esti¬ 


mate  of 


|Ys<*.y> 


cannot  be  made.  We  shall,  therefore,  use  as 


our  estimate  of  the  magnitude  of  Y_(x»y)  the  magnitude  of  the 

s 

estimate  ofYs(*»y)*  namely  |Yse(x»y)j* 


D.  NOISE  AND  SMOOTHING 

In  this  section  a  performance  index  which  is  an  average- 
signal-to-average-noise  ratio  is  defined,  and  expressions  for 
this  ratio  are  obtained  in  terms  of  useful  system  parameters. 

The  energy  received  from  a  unit  point  scatterer  is  given  by 

E(x,y)  3  J*|h(t;x,y)|  dt  (2-44) 

T 

which,  in  general,  is  a  function  of  x  and  y.  The  average  energy 
over  the  region  R  is  given  by 

E0  =  -•  JJ  J|h(t;x,y)j  dt  dx  dy  (2-45) 

R  T 


where  R  is  used  to  designate  both  the  region  occupied  by 
scatterers  and  its  area. 

It  follows  from  Eqs.  2-17  and  2-45  that 

oo 


EqR 


k=l 


(2-46) 


Furthermore,  if  SjJx  ,y;x'  ,y'  )  =  0,  it  follows  from  Eq.  2-JO' .that 

2  ,°° 

//|s(x-x«,y-y«>|  dx •  dy •  «  ^  |Pk(x,y) 

R  k-1 


2 
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for  x  and  y  inside  R.  Since  the  cross  section  of  S(x,y)  is  much 
smaller  than  the  area  R,  we  can  rewrite  the  above  expression  as 
follows: 

//|  S(x.y)  |  di  ^  pk! 

-  00  k=l 

2 

where  is  given  by  Eq .  2-43a.  Making  use  of  Eq.  2-29.  we 
obtain 

00  _ 

k  =  l 


It  is  convenient  to  define  the  following  sets  of  normalised 
quantities 


(2-48) 


£ 

n-1 


and 


k  oo 


(2-49) 


Xw 


n=l 


Applying  Eqs.  2-46  through  2-49  to  Eq.  2-43  «nd  redefining  terms 
gives 

N„  T 

(2-50) 


0  J 


R  EQ  <TR 


where 


(2-51) 


Because  O  is  much  smaller  than  R,  ¥  (x,y)  is  negligible  outside 

s 

R  except  in  a  region  that  (for  regions  R  of  reasonable  shape)  is 

negligible  relative  to  R.  Thus,  the  average  signal  intensity 


28 


s  fj\\i*-y'\  dx  dy  is  very  closely  approximated  by 
R 

°o  2 

S  ’  R  //|YsU,y)l  d*  dy  (2-52) 


-  00 


and  the  average-signal-to-average-noise  ratio  is 

_  _  S 

7}-  i c 


(2-53) 


Substituting  Eq.  2-27  into  Eq.  2-52  and  interchanging  orders 
of  integration  gives 
oo  oo 

S  “  ^  JJ  /TYR(x»,y«  )Ks#(x'-x",y»-y«)YR*(x".y")dx'dy,dxMy« 

'  (2-54) 


-  00  -  oo 


where 


oo 


Ksg(x*  -x"  ,y.'  -y"  )  *  JJ  S(x-x',y-y')  S* (x-x" ,y-y" )  dxdy  (2-54a) 


-  00 


YRU,y) 


(2 -54b) 


and  where 

Y(x,y)  for  (x,y)  inside  R 
0  otherwise 

Changing  variables  and  integrating  over  x*  and  y'  in  Eq.  2-54 
gives  „ 

S  "  R  Kss(u*v)  du  dv  (2-55) 


where 


-  oo 


00 


K\J/^(U»V)  ■  JJ  YR(x' ,y' )  YR*(x' -u,y-* -v)  dx‘dy'  (2-55a) 

-  oo 

whioh  is  the  autocorrelation  function  of  the  scatterer  density. 

Two  classes  of  scatterer  densities  for  which  explicit  results 
can  be  obtained  are: 
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1*  the  class  of  smooth  scattered  densities-- 

all  scatterer  densities  whose  autocorrelation 
function,  K^^(u,v),  does  not  change  apprecia¬ 
bly  over  a  region  of  sise  <r  centered  at  u  *  0 
and  v  ■  0 ; 

V 

2.  '  the  class  of  rough  scatterer  densities-- 

all  scatterer  densities  whose  autocorrelation 
function,  K^/(u,v),  in  the  vicinity  of  u  =  0 
and  v  =  0  is  substantially  confined  to  a  region 
much  smaller  than  ( T  .  (The  behavior  of  K^^u.v) 
for  u  and  v  outside  an  area  or  about  the  origin 
has  negligible  effect  on  S — see  Eq.  2-55«) 

If  we  note  that  K  (u,v)  is  a  spike-like  form  whose  cross 

ss 

section  is  approximately  or  ( see  Eq.  2-54a)  and  if  we  use  the  sub¬ 
scripts  s  and  r  to  refer  to  the  smooth  and  rough  class,  respec¬ 
tively,  then  Xq.  2-55  beoomes 

00 

Ss  *  Jl^fO.O)  JJ  Kss(u,v)  du  dv  (2-56) 

-  -  CO 

for  smooth  scatterer  densit ies  , . and 

oo 

Sr“RKss(0,0)  J/^u-v)  dudv  (2-57) 

-  00 

for  rough  scatterer  densities.  With  the  help  of  Eqs.  2-28,  2-29, 
2 -54a ,  2-55a,  and  2-54b,  Eqs.  2-56  and  2-5?  may  be  rewritten  as 
follows: 

S#  -  ^  //|^(x,y)|2  dx  dy  (2-58) 

R 

and 

sr  "  £-  |r  d*  dy|2  u-59) 

R 
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Substituting  Eqs.  2-50,  2-58,  and  2-59  into  Eq.  2-53  and 
redefining  terms,  we  obtain  for  the  smooth  and  rough  scatterer 
densities,  respectively , 


and 


.  _  !o  si 

7r  '  No  J 


(2-60) 


(2-61) 


where  IyI2  is  the  average  value  of  the  squared  scatterer  density 
1  —  12 

and  y  is  the  square  of  the  average  value  of  the  scatterer 
density;  that  is, 


|y|2  ■  £  JJ|Y(x,y)|2  dx  dy  (2-62) 

R 

and 

IyI2  =  |i  JJy(x« y>  d*  dy  J ^  (2-63) 

R 

In  order  to  provide  an  interpretation  for  J  in  Eqs.  2-60 
and  2-61 ,  suppose  ,  for  a  moment,  that  all  the  u-^'s  are  equal  for 
1  -  k  -  M  and  are  zero  otherwise.  The  sums  in  Eqs.  2-30,  2-48, 
2-49,  and  2-51  then  are  over,  at  most,  M  terms,  and  J  is  equal 
to  M--the  number  of  normal  modes  contained  in  h(t;x,y).  Next, 
suppose  that  the  values  of  the  1 s  are  approximately  constant 
for  1  -  k  -  M  and  decrease  very  rapidly  for  k  >  M.  If  has 
any  significant  value  for  k  >  M,  then  J  will  be  large  and  7)  and 
rjT  will  be  small  (see  Eqs.  2-60  and  2-6l).  As  a  result,  the  1  s 
must  be  confined  substantially  to  1  -  k  -  M.  and  then,  J  is  again 
approximately  equal  to  >1.  Thus,  in  general,  J  may  be  interpreted 
as  the  effective  numbqr  of  modes  that  may  be  utilised  without 
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undue  decrease  in  the  average-signal-to-average-noise  ratio. 
Despite  the  infinite  number  of  normal  modes  in  h(t;x,y),  the 
number  of  modes  that  are  usable  to  obtain  Y(x,y)  is  finite 

3  3 

because  the  high-order  modes  contain  little  energy  and  are  sub¬ 
merged  in  noise. 

If  the  smoothed  scatterer  density  is  to  be  an  exact  replica 
of  the  actual  scatterer  density,  then  (a)  S(x-x' ,y-y' )  must  be  a 
two-dimensional  impulse  6 (x-x1  ,y-y' ) ;  and  (b)  the  scatterer 
density  Y(x ,y)  must  contain  no  ambiguous  component  Yj_(x,y).  In 
order  for  condition  (b)  to  obtain  for  all  L f  unctions  Y(*  ,y)  . 
the  set  of  8k(x,y)  must  be  complete.  In  this  event  the  impulse 
may  be  regarded  as  the  limit  approached  by  the  sequence  of 
functions 


£ 
k  =  l 


.y' ) 


because  for  a  complete  set  of  9^(x,y) 


n — *  R  LEsf 


ek(x »y)  ®k*(x‘ ,y' ) 


dx'dy'  =  Y(x,y) 


The  corresponding  sequence  of  J's  may  be  evaluated  from  Eq.  2-51. 
giving  n 

n-i.2,  3.  ... 


k-1 


Since  the  Pk»  s  are  positive,  JR  >  or  Pn  > 

r  R 

finite  signal  energy  Bq,  the  series 


For  a 


'  oo 

XV 

k-l 
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converges  by  Eq.  2-46.  Hence,  according  to  Eq.  2-49,. 


n“l 


and  the  series 


oo 


must  also  oonverge.  Slnoe 


n=l 


diverges,  we  conclude  that  Jn  must  approach  infinity  with 
increasing  n.  Hence,  exact  reproduction  of  an  arbitrary 
scatterer  density ,y )  requires  an  infinite  number  of  usable 
modes . 

In  general,  as  the  form  of  S(x,y)  is  adjusted  to  shrink  cr 
toward  zero  in  an  effort  to  obtain  fine  resolution,  if  at  the 
same  time  ambiguities  in  the  determination  ofYg(x,y)  are 
avoided,  then  J  approaches  infinity  for  any  spike-like  shape 
S(x,y).  By  Eqs.  2-60  and  2-61,  the  average-signal-to-average- 
noise  ratios  T)s  and  7)  each  tend  to  zero  in  the  process.  We 
conclude  that  the  sharper  the  details  of  Y(x,y)  that  we  attempt 
to  estimate,  the  lower  the  average-signal-to-average-noise  ratios 
become.  As  a  result,  the  resolvable  cell  size  cr  must  be  chosen 
to  compromise  between  two  requirements: 

1.  (T  must  be  sufficiently  small  to  assure  that  the 
essential  details  of  Y(x,y)  are  not  smeared  out 
in  the  smoothing  process 

2.  (T  must  be  sufficiently  large  to  achieve  a  desired 
signal-to-noise  ratio 
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1.  THE  RECEIVER 

The  input  to  the  radar  receiver  is  r(t),  and  the  output 
at  the  end  of  the  observation  interval  T  is  a  two-dimensional 
display  of  ^a0(x,y)  (or  |^se(x»y)|)*  Frotn  E<ls*  2-33a,  2-34, 

2 -3 Ob,  2-40a,  and  2-40,  it  follows  that,  if 


r — '  Pv  ( x  ,  y ) 

Ws(x,yjt)  -  l.  i.m.  >  -a-— -  <Pk*(t)  (2-64) 

n"“*®k=l  k 


for  each  x  and  y,  then 

Jwa(x,y;t)  r  (t )  dt  “^se(*.y)  (2-63) 

T 

Thus,  the  reoeiver  is  described  by  the  linear  operation, 

f  Ws.(x,y;t )  r (t )  dt  • 

T 

if  \f/  (x  ,y )  is  the  desired  output.  As  long  as  J  is  not  infinite 
s  o 

and  cr  is  not  zero,  (see  Eq.  2-50)  is  finite  so  that,  as  a 
result  of  Eq.  2-43,  |n(x,y)|2  is  finite  almost  everywhere  in  R. 
Thus,  comparing  Eq.  2-64  with  Eq.  2-42,  we  see  that  the  right- 
hand  side  of  Eq.  2-64  converges  in  the  mean  as  a  function  of 
t ,  x ,  and  y . 

For  any  ®n(x»y)  that  is  nearly  constant  over  the  averaging 
region  of  the  smoothing  function  S(x,y),  Eq .  2-30a  yields 
pn(x,y)  «  6n(x,y).  Thus-,  if  the  0n(x,y),s  are  arranged  in  order 
of  decreasing  smoothness,  the  leading  terms  in  Eq .  2-64  will 
resemble  the  corresponding  terms  for  the  noiseless-case  receiver 
given  by  Eq.  2-25.  F or  a  given  l?l  (or  |^|  )  and  a  lower  bound 
on  I7r  (or  *)8)  »  Eq».  2-60  and  2-6l  indioate  that  the  lower  the 
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input  noise  level  NQ  or  the  greater  the  signal  energy  EQ ,  the 
closer  to  an  impulse  can  S(x,y)  be  made.  Thus,  as  the  input 
noise  level  decreases,  more  and  more  of  the  terms  in  Eq.  2-64 
will  resemble  terms  of  Eq .  2-25  for  the  noiseless-case  receiver. 

In  order  to  relate  the  above  results  to  the  conventional 

♦ 

theory  of  point -scatterer  detection  by  radar,  we  consider  a  point 
target  at  (xQ#yg)«  The  smoothed  scatterer  density  is,  by  Eq.  2-27, 

Ys(x,y)  =  S{x-x0,y-yQ) 

and  the  signal  intensity  on  the  radar  display  at  the  position  of 
the  target  is 


Ys(x0 


yP} 


|XZVX0’y0)  V^P^O* 

k 


where  we  have  made  use  of  Eq.  2-30.  With  the  help  of  Eq.  2-42, 
we  obtain  the  signal-to-noi se  ratio  at  the  position  of  the  target 


^s^o-yp^l 
n<*0»y0)  1 2 


1)  .kv^o-vlf 


PkUo.yp) 


Applying  the  Schwarts  inequality  to  the  numerator  gives 

Y3(*o.y0)l2 

n(x0,y0)|2 

where  we  have  also  made  use  of  Eqs.  2-44  and  2-17.  The  maximum 
signal~to-noise  ratio  is  achieved  when 


»o  XI 


®k <*o ' 


y0J 


N. 


Pk(xo,yp) 


c(xp*yp)  ek(xp,yp) 


(2  -66) 


where  c(x0,yQ)  is  an  arbitrary  function  of  xQ  and  y0  and  is 
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independent  of  k.  Thus,  aooording  to  Eqs.  2-64  and  2-66,  the 
reoeiver  which  maximises 

1  ’Vl2 

n<*o’yo)|2 


for  all  x«  and  yft  inside  the  region  R  is  given  by 


n 

W0(x,y;t)  -  c(x,y)  X.i.m.  \  M-k  6k(x,y)  <p  *(t)  (2 -67) 

where  WQ(x,y;t)  is  used  to  designate  the  Ws(x,y;t)  obtained  when 
Pktx.y)  is  given  by  Eq.  2-66.  The  resulting  smoothing  function, 

which  may  not  have  a  spike-like  form  5<x-x' ,y-y' )  is  given  by 

oo 


S(x,y;x' ,y> )  ■  c(x 


^k  0k(x*y)  V(x’’y,)  (2_68) 


k=l 


Comparing  Eqs.  2-67  and  2-68  with  Eq .  2-17,  we  obtain. 

Vq(x  ,y  ;t )  *»  c(x,y)  h*(t;x,y) 

and 

S(x,y;x' ,y')  *  c(x,y)  Jh(t;x,y)  h*(t;x',y' )  dt 

T 


Thus,  aside  from  the  multiplying  function  c(x,y),  the  receiver 
reduces  to  the  matched  filter;  the  smoothing  function  reduces  to 
a  generalized  form  of  Woodward's  uncertainty  function;  and  the 
output  signal-to-noise  power  ratio  reduces  to  the  energy  ratio 
B(x,y)/NQ  of  conventional  detection  theory. 


CHAPTER  III 


DELAY-DOPPLER  RADAR  THEORY 

A.  THE  SC ATTERER -DENSITY  MODEL 

It  is  assumed  in  this  chapter  that  x  is  time  delay  and  y 
is  Doppler  angular  frequency.  It  is  further  assumed  that,  if 
f(t)  is  the  complex-video  representation  of  the  transmitted  wave¬ 
form,  the  radar  return  from  a  point  scatterer  at  delay  x  and 
Doppler  angular  frequency  y  may  be  expressed  as 

h(t;x,y)  =  f(t-x)  e~^ty  (3“D 

As  a  result  of  the  superposition  of  electromagnetic  fields 
in  linear  media,  the  radar  return  due  to  any  transmitted  wave 
f(t)  can  always  be  represented  by 

oo 

s(t)  =  J g(x,t)  f(t-x)  dx  (3-2) 

-  00 

I 

where  s(t)  is  the  received  waveform  and  g(x,t)  is  independent 
of  f(t)  and  depends  only  on  the  antenna,  the  scatterer  distribu¬ 
tion,  and  their  relative  positions  and  motion. 

On  physical  grounds  it  may  be  assumed  that  g(x,t)  is  an 
integrable-square  function  of  time  over  any  finite  time  interval 
L  so  that  the  integral 

Y(x.y)  -  ^  Jg(x,t)  ejty  dt  "  (3-3) 

L 


3  6 
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exists.  Hanot ,  by  Fourier  transform  theory,  we  have 

oo 

g(x,t)  -  J  e“^ty  dy  (3-JO 

-  oo 


for  t  inside  the  interval  L.  Substituting  Eq.  3 -k  into  Eq.  3-2 

00 

s(t)  -  JJ  \fr(x,y)  tit-x)  e'jty  dx  dy  (3-5) 


-  00 

for  t  inside  L.  The  interval  L  is  assumed  to  be  sufficiently 
long  so  that  any  observation  interval  T  of  interest  to  the  radar 
user  falls  inside  L  and  Eq.  3~5  is,  therefore,  valid  for  t  inside 
T.  Equation  3-5  is  the  same  as  Eq.  2-3  if  h(t;x,y)  is  defined  by 
Eq.  3-1 .  Clearly  h(t;x,y)  depends  only  on  the  radar;  and  by  Eq. 
3-3.  \f/lx,y)  ,  like  g(x,t),  depends  only  on  the  scatterers  and  the 
radar  antenna.  To  eliminate  dependence  of  <//(x,y)  on  the  radar, 
we  assume  a,  priori  that  the  scatterers  are  confined  to  a  region 
of  uniform,  unchanging  antenna  illumination.  Thus,  the  condi¬ 
tions  required  of  the  coordinates  x  and  y  in  Chap.  II,  Sect.  A, 
are  satisfied  and  determination  of  ^(x,y)  can  be  considered  the 
radar  objective  in  a  study  of  performance  obtainable  with  differ¬ 
ent  transmitted  waveforms  and  receiver  processing  functions. 

Note  that,  in  view  of  the  choice  of  h(t;x,y)  (see  Eq.  3-1)  ,  it 
is  to  be  expected  that,  subject  to  the  usual  restrictions  (such 
as  wide  antenna  beamwidth,  narrow  transmitter  bandwidth,  negli¬ 
gible  target  acceleration),  |^(x,y)|  may  be  regarded  as  an 
approximate  map  of  the  distribution  of  objeots  throughout  the 
x-y  plane. 
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B.  PROPERTIES  OF  TRANSMITTED  WAVEFORMS 

It  is  assumed  in  the  remainder  of  this  chapter  that  Y(x,y) 
is  known  a  priori  to  be  zero  for  x  outside  an  interval  D,  extend¬ 
ing  from  x  =  x^  to  x  =  Xg.  It  will  be  shown  that,  in  consequence, 
certain  classes  of  waveforms  need  not  be  considered  in  evaluat¬ 
ing  radar-performance  possibilities. 

Consider,  f or  a  moment,  a  nonperiodic  transmitted  waveform 
f(t)  in  relation  to  the  time  scales  for  f(t)  and  s(t)  in  Fig.  2. 
Since  x^  is  the  minimum  delay  and  x ^  is  the  maximum  delay,  it 
follows  from  Fig.  2  that  any  transmission  outside  the  interval 
T^  will  make  no  contribution  to  s(t)  inside  the  observation 
interval  T  (which  for  convenience  is  centered  at  t  =  0  ) .  Thus, 
the  transmitted  waveform  could  just  as  well  have  been  a  periodic 
waveform  of  period  T^.  It  is  therefore  concluded  that  only 
periodic  transmitted  waveforms  need  be  considered. 

Next,  suppose  that  a  periodic  waveform  of  period  P,  which 
is  less  than  the  length  of  the  interval  D,  is  used  as  the  trans¬ 
mitted  waveform.  Since  f(t)  =  f(t  +  P),  it  is  clear  from  Eq.  3-5 
that  any  change  inY(x,y)  coupled  with  the  negative  of  that 
change  shifted  by  P  along  the  x-axis  ,  and  such  that  both  changes 
are  confined  to  D,will  not  affect  the  radar  return.  Hence,  if 
P  <  D,  an  ambiguous  component  of  the  scatterer  density  will 
arise.  Waveforms  of  period  less  than  D  are,  therefore,  of  no 
interest  if  unambiguous  determination  ofY  (x,y)  is  desired.  In 
the  following  section  we  consider  waveforms  of  period  P  =  D. 
Attempts  to  treat  waveforms  of  period  P  >  D  have  been  unsuccess¬ 
ful  --it  has  not  been  possible  to  determine  whether  or  not  limita-  . 
tions  present  for  P  ■  D  can  be  alleviated  for  P>D. 


I\>|H 
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C.  AMBIGUITIES  OF  RANGE-DOPPLER  RADAR  FOR  P  =  D 

In  this  section  we  consider  transmitted  waveforms  of  period 
P  equal  to  the  length  of  the  scatterer-occupied  interval  D.  We 
shall  be  interested  in  establishing  conditions  under  which  a 
sufficiently  complete  set  of  6^(x,y)'s  can  be  obtained. 

The  class  of  waveforms  under  consideration  can  be  represent 
ed  by* 

fit)  -Y2**  •J2’,kt/D  (3-6) 

k 

where  the  fk  are  the  Fourier  coefficients,  and  D  is  the  length 
of  the  interval  D.  Furthermore,  for  t  inside  the  observation 
interval,  we  have 

h(t ;x ,y )  “^h^x.y)  *^2lTlct/T  (3-7) 

where 

hk(x,y)  =  —  Jh ( t ; x , y )  e"J2TTkt/T  dt  (3-8) 

T 

Also  note  that,  without  destroying  any  useful  information  in  the 
radar  return,  the  observation  time  T  can  always  be  extended  so 
that 

T 

jj-  =  I  =  Positive  integer  (3 -9) 

From  Eqs.  3-1,  3-6,  3 -7,  3-8,  and  3-9 ,  it  follows  that 

h(t;x,y)  e-J2’,tWD 

n,£ 

Equation  3-5  and  3-10  indicate  that  any  component  of  the 
scatterer  density  which  is  orthogonal  to  for 

♦Henceforth,  it  is  understood  that  all  summation  indices  run 
from  -00  to  +00,  unless  otherwise  indicated. 


sinU(IA-n),  -  yT/2  J  J2rr«it/T  /o_10) 
iT(If-n)  -  yT/2  ®  10' 
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k  will  make  no  contribution  to  the  radar  returns.  Hence,  the 
component  that  does  contribute  to  the  radar  return  must  be 
expressible  as 

YT(x,y)  •  ^rVT(x,2TTk/T)  3  (3-11) 

k 

It  has  already  been  assumed  that  Y(x,y)  =  0  outside  an  x 
interval  of  length  D.  It  is  now  assumed  also  that  the  samples 
YT(x,2Ttk/T)  are  known  a_  priori  to  be  zero  for  k  <  k^  and  k  >  kg. 
For  T  >>  D  this  assumption  is  often  nearly  equivalent  to  assuming 
Y(x,y)  =  9  outside  a  band  of  Doppler  angular  frequencies  y, 
extending  from  y^^  a  2nk1/T  to  y2  “  2frk2/T,  The  width  of  this 
band  is  given  by 

B  =  y2  -  yx  =  2Ti(k2-k1)/T  (3-12) 

Equation  3~H  can,  therefore,  be  rewritten  as  follows: 

(3-13) 

Henceforth,  any  scatterer  density,  whose  Yq^x.y)  component  can  be 
expressed  as  a  finite  series  in  the  form  of  Eq.  3“13  with  B 
determined  by  Eq.  3-12  and  which  is  zero  for  x  outside  D,  will 
be  referred  to  as  a  BD  scatterer  density. 

Observe  that,  in  the  absence  of  any  additional  a  priori 
inf  ormation  ,  YT(  x  ,y )  is  the  only  component  of  the  scatterer 
density  that  can  possibly  be  deduced  from  the  radar  return.  We , 
therefore,  wish  to  find  the  maximum  value  of  B,  if  any,  for 
which,  in  the  absence  of  noise , Yt^x *y)  can  be  deduced  from  the 
radar  return.  We  also  wish  to  determine  what  restrictions  on 
transmitted  waveforms  are  needed  to  permit  deduoing  aYT(x,y) 


of  maximum  B 
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Xt  is  convenient  to  define  a  set  of  coefficients 

when  k/l  is  an  integer 


Si 


otherwise 

and  to  rewrite  Bq.  3-10  in  the  form 

h(t  ;x  ,y)  ■  h0(t;x,y)  +  h^tjx.y) 

h0(t ;x,y) 


(3-14) 


(3-15) 


where 


n  m=k. 


-  j2fr(n+m)x/lD  sin  (rrm  -  vT/2 )  j2rrnt/T 
gn+m  9  tto  -  yT/2 

(3-16) 


and  h^(t;x,y)  has  the  same  form,  but  the  index  m  runs  over  values 
outside  the  (^^2)  interval.  It  follows  from  Bq.  3-13  that 

00 

jYT(x,y)  hjUjx.y)  dy  -  0  (3-1?) 

-  00 

Thus,  h^tjx.y)  is  of  no  concern  to  us  since,  for  a  BD  scatterer 
density,  it  makes  no  contribution  to  the  radar  return. 

Equation  3-16  may  be  written 


where 


hQ(t  ;x,y) 


21  “n 


(3-18) 


n+k2 

^  N2  °'l9> 

i*n-*-k1 


8  (x  y)  -  V"  g*  j2n(m+n)x/lD  HpW m-  yT/2)  <3_20) 

n'*,y/  \i  /  .  6m+n  "  trm  -r  yT/2  ' 

m=ki 

and 

.It )  =  -i_  .J2"n*/T  (3-21) 

Note  that,  if  any  M-n  *  0  when  calculated  by  Bq.  3-19,  then  the 


corresponding  term  is  absent  from  series,  Eq.  3-18,  so  that 
there  is  no  need  to  consider  indeterminates  in  Xq.  3-20.  Vith 
the  help  of  Eq*  3-14,  it  oan  be  ehown  that 

n 

(3-22) 

n 

where  the  region  R  is  bounded  by  the  interval  D  along  the  x-aiis 
and  extends  from  -  co  to  +  oo  along  the  y-axis.  Thus,  the  0n(x,y) 
form  an  orthonormal  set  over  R.  Furthermore,  the  <pn(t)  form  a 
complete  orthonormal  set  over  T. 

By  Eqs.  3‘1  and  3-5.  the  discussion' above  Eq.  3“H.  Eq.  3-1?. 
and  becauseYT(x  ,y )  *  0  for  x  outside  D,  we  have 

oo 

s(t)  «  J  JVT(x.y)  h0(t;x,y)  dx  dy  (3-23) 

D  -  oo 

for  t  inside  T.  Equations  3-23  and  3-18  are  analogous  to  the 
basic  equations  of  the  general  theory  (Eqs.  2-3  and  2-17).  The 
only  change  is  that  the  f unct ion YT (x ,y ) ,  instead  of  being 
identically  zero  outside  an  interval  B  (to  yield  a  finite  region 
R) ,  is  a  sin  x/x  interpolation  of  samples  spaced  at  intervals 
2tt/t  along  y;  and  the  sample  values  are  zero  outside  B.  The 
properties  of  orthonormal  expansions  obtain, as  in  the  general 
theory  (see  Eqs.  2-22  through  2-23f)  ;  theref  ore ,  YT(  x  ,y)  is  the 
sum  of  a  component  Y||(x,y)  expressible  as  a  linear  combination 
of  the  6n(x,y)  and  an  ambiguous  component  YjJx.y)  that  makes  no 
contribution  to  s(t)  and,  hence,  cannot  be  determined  by  the 
radar.  The  necessary  and  sufficient  condition  for  unambiguous 
determination  ofYT(x,y)  is  therefore  thatYT(x,y)  is  known 
a  priori  to  be  a  linear  combination  of  the  6n(x,y).  Ve , 


if  m  = 


if  m  j 
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therefore,  wish  to  determine  the  maximum  value  of  Q  for  which 
theN^lx  ,y)  component  of  any  BD  scatterer  density  can  be  expand¬ 
ed  in  terms  of  fl  (x,y). 

It  is  convenient  to  express  each  of  the  samples  Y(x ,2nk/T) 
in  Eq..  3-13  as  a  Fourier  series  over  the  interval  D.  Hence, 

k2 

YT(x,y)  -  ^  y>  (3-2U) 

m  n=k^ 

where 


u  (x  v)  =  eJ2lTmx/D  sin  (ntr  -  yT/2) 
0  rm  -  yT/2 


mn 


(3-25) 


and  where  the  Y  are  the  appropriate  Fourier  coefficients.  In 
mn 


order  to  be  able  to  expand  H/^,( x , y)  as  a  linear  combination  of 
the  6  (x.y).  it  is  necessary  and  sufficient  that  each  of  the 
orthogonal  functions  umrj(x,y)  be  representable  as  a  linear  super¬ 
position  of  the  8v(x,y);  that  is,  it  is  necessary  and  sufficient 


that 


u»nU-’r>  ‘  Y2* 


k‘“n)  ek(x.y) 


(3-26) 


where 


Uk(mn)  =  //umn(x’y)  V(x’y)  dx  dy  (3 -26a ) 

R 

for  all  m  and  for  k^  -  n  -  k2« 

Substituting  Eq.  3“26a  into  Eq.  3-26  and  making  use  of  Eqs. 
3-25,  3-20,  and  3-l4,  we  find  that  it  is  necessary  and  suffi¬ 
cient  that 


g 


2ttD 


Im 


Im-n  t =k. 


E  gIm-n+f 


e 


j2tt  (Im+i-n  )x/lD  sin  -  yT/2 ) 

tt l  -  yT/2 

(3-27) 


Thus,  if  Eq.  3-26  is  to  be  a  oorreot  representation  for  umn(* »y ) 
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as  Riven  by  Eq.  3-25,  it  is  neeessary  and  sufficient  that 


2ttD 

4-„  *In  gi™-n  I0  -  ■ 

for  each  m  and  for  k^  -  n  k2 .  Implied  in  Eq.  3-28  is  that  for 
no  k  can  ^  ■  0.  Observe  that,  if  any  ^  8  0,  then  the  corres¬ 
ponding  term  would  be  missing  from  Eq.  3-26  and  the  right-hand 
member  of  Eq.  3-27  would  be  zero  for  the  corresponding  value  of 
Im-n  so  that  Eq.  3-27  could  not  be  satisfied. 

The  above  necessary  and  sufficient  conditions  oan  be 
restated  as  follows: 


-I1 

Im+A-n 


when  t  -  n 


(3-28) 


Vn 
V  2ttD 


^  0  for  eaoh  m  and  k  *  k2  (3-29) 


«Im«im+t-n  *  0  for  «<»>"•  'cl'‘'k2'  kl  *  "  *  k2  .  <  #  " 


(3-30) 


If  we  let  k  *  l  -  n,  Eq.  3-28  becomes 

gIm  eIm+k  *  °  f°r  each  “•  -<k2  -ki)  (k2  ”  kl }  .  k  /  0 

(3-3D 

Equation  3-29  requires  that,  for  all  m,  g^m  /  0.  Therefore, 
referring  to  Eq.  3-31,  we  conclude  that,  in  order  to  satisfy 


Eqs.  3-29  and  3-31.  it  is  necessary  that 

k2  -  k2  <  I 


(3-32) 


which,  with  the  help  of  Eqs.  3-9  and  3“12,  Eq.  3"32  may  be 


restated  as  follows: 


BD  <  2tt 


(3-33) 


The  maximum  value  of  k2~k1,  which  satisfies  Eq.  3-32  is 


given  by 


k2  -  -  X  -  1 


(3-31*) 


4  6 


which  value  of  kg  -  k  represents  the  worst  possibility--if  an 
^adequate  set  of  6n(x,y)  can  be  achieved  for  kg  -  ^  =  I  -  1 ,  then 
it  oan  also  be  achieved  if  fewer  samples  x  ,2TTk/T)  are  permit¬ 
ted  by  ^  priori  knowledge  to  differ  from  zero.  Equation  3~34 
implies  that,  for  any  integer  n,  a  unique  integer  p  can  be  found 
such  that 

n  +  kl  *  Tp  "  n  +  k2  (3-35) 


Applying  Eq.  3-35  to  Eq.  3-19  and  making  use  of  Eq.  3-14,  we 
find  that 

2 


=  2ttD 


for 


Vkz'"'Vkl  <3-3«) 


In  other  words  for  every  transmitted  harmonic  f  ,  there  are 
kg  -  k^  +  1  =  I  equal  characteristic  values  ;  and  since  no 
fp  =  °*  then,  for  every  n,  n-n  j-  0.  From  Eqs.  3-36  and  3-14,  it 
follows  that  Eq.  3-29  is  satisfied.  Also,  as  a  result  of  Eq. 
3-34,  Eq.  3_30  is  satisfied.  Thus,  we  are  assured  that,  if 
(a)  for  all  p,  f  /  0;  and  (b)  Eq.  3-33  holds,  the  set  of  G^x.y) 
is  sufficient  to  expand  the  Y^(x,y)  component  of  a  BD  scatterer 
density,  and  unambiguous  determination  of Y^(x.y)  is  possible. 

The  set  of  8n(x,y)  is  obtained  by  using  Eq.  3-35,  together  with 
Eq.  3-14  in  Eq.  3-20 .  Hence, 


en(x,y) 


/ 


T 

2ttD 


j2TTpx/D  sin  [tt(  Ip_n)  -  vT/2. 

^(Ip-n)  -yT/2 


(3-37) 


for 


I  -  k_  -  n  I  -  k. 


so  that,  for  every  transmitted  harmonic  f  ,  there  are 
kg  ”  k^  +  1  =  I  characteristic  functions  &n(x,y). 

For  a  transmitted  wave  of  period  P  =  D,  we  have,  therefore, 
shown  that,  in  order  to  obtain  unambiguous  determination  of 
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^T(x,y)  for  a  BD  scatterer  density,  It  is  necessary  and  suffi¬ 
cient  that 

(a)  the  transmitted  ware  oontains  all  its  harmonies 
(no  fk  ■  0) 

(b)  kg  -  <  I,  or  BD  <  2tt  --  that  is,  the  scatterer 

density  must  be  such  that  the  samples  ^(x, 2nk/T) 
are  zero  for  k  <  k^  and  k  >  kg 

It  is  interesting  to  note  that  condition  b  corresponds  to 
the  well-known  limitation  on  maximum  unambiguous  delay  and  maxi¬ 
mum  unambiguous  Doppler  frequency  encountered  in  periodic  pulse 
radar  systems. 


D.  NOISE  AND  SMOOTHING 

It  is  assumed  in  this  section  that  Eq.  3-33  is  satisfied, 
that  the  period  of  the  transmitted  waveform  is  equal  to  D,  and 
that  all  harmonics  are  present.  Thus,  we  are  assured  that  there 
is  no  ambiguous  component  of  \jr^(x,y).  In  analogy  to  Sect.  D  of 
Chapter  II,  a  performance  index  which  is  an  average-signal-to- 
average-noise  ratio  will  be  defined.  Expressions  for  this  ratio 
will  be  obtained  in  terms  of  useful  system  parameters. 

Suppose  that  a  smoothing  operator  of  the  following  form  is 
chosen 

S( x ,y )  =  a (x )  b ( y )  (3-38) 

where  both  a(x)  and  b(y)  could  be  some  pulse-like  forms  centered 
at  x  =  0  and  y  =  0 ,  respectively.  The  width  of  a(x)  is  much 
smaller  than  D,  and  the  width  of  b(y)  is  much  smaller  than  B. 

We  oan ,  therefore ,  write 
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.(.)  (3-39) 

k 

In  order  to  assure  that  S(x-x',y-y')  is  expressible  in  terms  of 
the  @n(x,y)  (see  Chap.  II,  Seot.  C)  ,  it  is  necessary  to  assume 
that  the  Fourier  transform  of  b(y)  is  eero  outside  an  Interval 
T.  Henoe, 

27  J  ^<y-u)  J~uT'/2^  du  =  b<y)  (3-4o) 

-  00 

To  apply  the  theory  of  Chap.  II,  Sect.  C,  we  note  that,  in 
analogy  to  Eq.  Z-JOa  ,  we  have 


VAS 

?n  ( *  •  y )  =  f  J  S(x-x',y-y')  6n(x' ,y' )  dx'dy'  (3-41) 


D  -  oo 


Thus,  using  Eqs.  3-38,  3-39,  and  3-40  in  Eq.  3-41  gives 

* 

P„U.y)  -yS  r^r  *p  b[y  -  2*<Ip-n)] 


for  Ip-lc^n^Ip-k, 


(3-42) 


Furthermore,  by  analogy  to  Eq.  2-42  and  with  the  help  of  Eqs, 


3-40,  3*42,  and  3 “36,  we  obtain 


In (x  ,y  )|  2  s  Y  7^  y  [  |b(y  -  2mn/T)|  (3-43) 


,2  k2 


The  assumption  resulting  in  Eq.  3-40  also  yields 


b (y-u )  b(y  -  2TTm/T) 


sin  ( Tim  -  uT/2 ! 
mn  -  uT/2 


(3-44) 


oo  _ _ ^ 

J  | b (y  — u )  | 2  du  ■  y  y  jb(y  -  2TTm/T)j 
-oo  m 


so  that 


(3-45) 
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Comparing  this  expression  with  Eq.  3-43 ,  we  obtain 

|n(x.y)|2  *  Kr 

where 

N  .  a  2  oo  2 

nr  "  57  XI  f2  /lb<y,l  dy  (3‘46) 

p  P  -  00 

Since  b ( y )  is  a  narrow  spike,  much  narrower  than  B,  we  see  that, 
for  y  inside  B,  and  not  too  near  the  edge  of  B,  we  have  approxi¬ 
mately  , 

|n(x,y)|2  -  for  y  inside  B  (3“47) 

Furthermore,  for  y  appreciably  outside  B,  jn(x,y)|2  is  nearly 
aero. 

The  quantity  N^,  given  by  Eq.  3“46,  cannot  be  defined  the 
same  way  as  in  the  general  theory  (see  Eq.  2-43)  because  the 
region  R  in  this  chapter  is  infinite.  It  is  interesting,  though, 
that 

oo  _ 

bd  f  J  |n(x *yV| 2  dx  dy 

D  -  oo 

leads  exactly  to  above;  and  to  that  extent,  BD  is  the  "effective" 

area  of  R  and  is  the  "effective"  average  of  the  mean-square  noise. 

Because  the  region  R  considered  here  is  infinite,  the 'average 

signal  intensity  also  cannot  be  defined  as  in  Chap.  II  (see  Eq.  2-52). 

Nonetheless,  the  average-signal-to-average-noise  ratio  calculated 

from  Eqs.  2-42,  2-52,  and  2-53, 

00  2 
//|YsU.y)|  dxdy 

77“  . . .  (3-48) 

jy|n(*.y)j2  dx  dy 
R 


i 
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exists  inasmuch  as  the  factor  l/R  is  cancelled  in  the  ratio. 

By  a  procedure  similar  to  the  one  used  in  obtaining  Eq.  2-57,  j 

we  find  that,  for  a  rough  scatterer  density 
00  00 

//  |^s(x’y>|2  dxdy  "  KSS(0*0)  JJKWU,v)  dUdV 

-  00  -  oo 

whioh,  with  the  help  of  Eqs.  2-54a  and  2-55a  ,  becomes 

co  oo 

ff\V  x , y ) J  2  dxdy  =  //I-  ,y)|2  dxdy  |//t(x.y)  dxdy 

-  oo  -  oo  R  (3 -49) 

A  similar  relation  can  be  obtained  for  the  smooth  case;  but  we 
shall  confine  our  attention  to  the  rough  case,  which  appears  to 
be  of  greater  physical  significance.  It  follows  from  Sqs.  3-43, 

3-9.  3-38,  and  3-49  that  Eq.  3-48  can  be  written  as 

oo 

/|a(x)|2  dx 

V-  - j  \fj tU.y)  dxdy |  (3-JO) 

No  zl  k/frl  R 

P 

It  is  interesting  to  define  the  number  of  useful  modes,  a 
quantity  analogous  to  J  given  by  Eq. ,2-51  in  the  general  theory. 

The  normalized  quantities  ci^.  and  pk  can  be  defined  as  in  Eqs. 

2-48  and  2-49,  except  that  pk2  in  Eq.  2-48  cannot  be  defined  for 
an  infinite  region.  However,  only  the  ratio  of 


i 

f 

l 

t 

) 


is  needed;  henoe ,  we  define 


//|pk(x.y)|2  dxdy 


R 


y~  //i,n<*.H2  “i,iy 


n  R 


Making  use  of  Eqs.  2-49,  2*51,  3-36,  3-39.  3-42  and  the  above 
definition,  va  obtain 


<3-51) 


whora  PQ  is  tha  average  transmitted  power  given  by 

vEW  <3'52 

n 

Note  that,  as  a  result  of  Eqs.  3-1.  3-6,  3-9.  and  3-52,  we 
can  write 

BD  ///h<t;X,y)|2  dt  d*  dy  *  P0  T 
B  D  T 


regardless  of  the  length  of  the  interval  B.  Comparing  the  above 
result  with  Eq.  2-45,  we  have 

E0  =  P0  T  (3-53) 

With  the  help  of  Eqs.  3-51,  3-53,  and  3-39,  we  oan  rewrite  Eq.  3-50 
as  follows: 

7lm  jj-^7  |//Vu.y)  dx  dy|  (3-54) 

0  R 


whioh,  if  we  take  into  account  Eq.  2-63,  agrees  with  Eq.  2-61. 
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E.  SYNTHESIS  OF  THE  TRANSMITTED  WAVEFORM 

In  this  section  we  obtain  transmitted -waveform  requirements 
which,  for  a  given  transmitted  average  power  and  the  smoothing 
function  specified  by  Eq.  3-38,  yield  the  maximum  output  signal- 
to-noise  ratio. 

As  long  as  the  f  1 s  are  different  from  zero,  the  smoothed 

scatterer  density  \{/  (x,y)  is  independent  of  the  f  1  s  so  that  we 

®  _  P 

can  choose  a  set  of  f  1 s  to  minimize  |n(x,y)|2,  as  given  by 
Eq.  3“**3  without  affecting  the  signal  intensity.  Thus,  the  set 


of  f  '  s  which  minimizes  |n(x,y)|  also  maximizes  both  the  average- 
signal-to-average-noise  ratio  as  well  as  the  local  signal-to-noise 


ratio  |\^s(x,y)j  2/|n(x,y)|2  . 


It  is  clear  from  Eq.  3“^3  that  the  larger  the  Jf^j  1  s  the 
smaller  n (x  ,y ) I  2  .  However,  since  the  average  transmitted  power 


PQ  must  remain  finite,  we  need  to  impose  Eq.  3-52  as  a  constraint. 
Applying  the  method  of  Lagrange's  multipliers,  we  minimize  Eq.  3~*+3 
under  the  constraint  imposed  by  Eq.  3-52.  The  result  is 


(3-55) 


Thus ,  once  the  smoothing  operator  has  been  selected,  the  optimum 
distribution  of  transmitted  power  throughout  the  spectrum  is 
determined  by  Eq.  3_55«  The  phase  angles  of  the  fn's  are  arbi¬ 
trary  so  that  it  might  be  possible,  in  principle,  to  minimise 
the  peak  power  requirements  by  adjusting  the  phases. 

If  a(x)  is  a  realizable  autocorrelation  function,  and  if 


the  optimum  transmitted  waveform  is  used,  then  Eq.  3 -51  yields 
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■afi&L 


00 


<3-5*> 


y  |*  (x )  |  2  dx 


-  oo 

where  we  have  also  made  use  of  Eqs.  3“ 39  and  3-55  •  The  factor 
.oo 


/QU  -  - 

|a(x)|Z  dx/az(0)  is  a  measure  of  the  width  of  a(x)  or  a 

-  oo 

measure  of  the  reciprocal  bandwidth  of  the  spectrum  of  a(x). 

However,  by  Eq.  3-55,  a(x)  is  the  autocorrelation  function  of 

2  f 00  2 

the  transmitted  wave  so  that  a  (0)/  I  |a(x)|  dx  is  a  (somewhat 

-  co 

unusual)  measure  of  the  bandwidth  of  the  transmitted  wave;  and, 
in  this  special  case,  J  reduces  to  a  radar  time -bandwidth  product. 

The  average-signal-to-average-noise  ratio  for  the  optimum 
transmitted  waveform  is  given  by  Eq.  3~5^  with  J  evaluated  as 
the  time-bandwidth  product  of  Sq.  3 -56. 


W.  THE  PROCESSING  OF  THE  RADAR  RETURN 


The  receiver  weighting  function  corresponding  to  the  choice 
of  smoothing  function  given  by  Eq.  3-38  is  obtained  in  the  same 
way  as  in  the  general  theory  (Chap.  II,  Sect.  F).  Equation  2-6l 
results,  except  that  the  summation  index  now  runs  from  -oo  to 
+  oo.  Substituting  Eqs.  3“21,  3-3$,  and  3-41  into  Eq.  2-64,  we 
obtain 


I  -k. 
J2_  1 


(3-57) 


n-Ip-k2 


Substituting  m  for  Ip- g,  making  use  of  Sq.  3-9,  and  rearranging 
terms  gives 
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m=k1  _ 

(y-2™/T>  .J2"*‘/Tj 


(3-58) 


It  is  convenient  to  define  the  Courier  transform  of  b(y) 


B(t) 


-  /  b(u) 


(3-59) 


from  which  it  follows  that 


w 

B(t)  e^y"  J  b(y-v)  e^vt  dv 


(3-60) 


Making  use  of  Eq.  3-44  in  Eq.  3-60,  we  obtain 

B(t)  .•W-f^blyin./ll  .J2™*/T  (3-61) 

m 

for  t  inside  the  interval  T,  If  we  recall  that  the  width  of  the 
spike-like  form  b{y)  is  much  narrower  than  B,  we  see  that,  for  y 
inside  B,  and  not  too  near  the  edge  of  B,  Eq.  3~6l  may  be  approxi¬ 
mately  rewritten  as  follows 

k2 

Bit)  .-•»**  -  &Y2  b(y-z™/x»  .J2™*/1  (3-62) 

ra=k^ 

Using  Eq.  3-62  in  Eq.  3-58  gives 

Ws,x.y:t)  .  .-■l2T’P<t-*>/D  Bit)  ,‘ir*  (3-63) 

P  P 

Suppose  that  the  spike-like  form  a(x)  is  a  realisable  auto¬ 
correlation  function  so  that  all  the  an  are  positive.  Then,  if 
the  optimum  transmitted  waveform,  specified  by  Eq.  3 -55.  is  used 


in  Eq.  3-63 


*9(x,y;t)  s|^B(t)  f*(t-x)  ejty 


(3-64) 
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where  we  have  also  made  use  of  Eq.  3 "39  and  Eq.  3-6.  If  B(t) 
is  constant  throughout  T,  so  that  b(y)  specifies  only  the 
unavoidable  smoothing  associated  with  a  finite  observation 
interval,  the  weighting  function  in  Eq.  3-60  is  the  matched 
filter  of  conventional  radar  theory. 


ft.  APPLICATION  TO  THE  MAPPING  OF  THE  LUNAR  SURFACE 

An  interesting  application  of  radar  to  dense-scatterer 
distributions  is  encountered  in  the  problem  of  radar  mapping  of 
the  lunar  surface  in  range  and  range-rate  (or  delay  and  Doppler- 
frequency)  coordinates.  Such  a  mapping  has  been  carried  out  by 
the  Millstone  Radar.11 

Of  importance  to  this  mapping  is  the  libration  of  the  moon  — 
an  apparent  angular  vibration  of  the  moon  as  viewed  by  an  observ¬ 
er  on  the  earth.  The  major  component  of  libration  arises  in  the 
following  manner:  The  same  hemisphere  of  the  moon  is  continually 
pointing  toward  the  center  of  the  earth.  Thus,  an  observer  on 
the  surface  of  the  earth  sees  slightly  different  portions  of  the 
moon  as  the  earth  rotates  to  place  the  moon  in  different  posi¬ 
tions  above  the  observer's  horizon.  As  a  result  of  this  effect, 
the  moon  appears  to  rotate  on  its  own  axis  through  a  small  angle 
when  viewed  from  a  point  on  the  earth's  surface.  Different 
points  on  the  surface  of  the  moon  will,  therefore,  have  differ¬ 
ent  range  rates.  Thus,  a  time-delay  and  a  Doppler-f requency 
coordinate,  together  with  a  scatterer  density  \j/(x,y),  may  bo 
assigned  to  eaoh  point  on  the  surfaoe  of  the  moon. 
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In  the  Millstone  experiment  an  estimate  of[i//(x,y)|  ,  where 

x  is  time  delay  and  y  is  Doppler  frequency,  was  obtained.  The 
estimate  may  be  assumed  to  differ  from  a  physical  map  of  the  lunar 
surface  for  reasons  pointed  out  in  Chap.  II,  Sect.  A.  Furthermore! 
aside  from  equipment  limitations,  noise  and  smoothing  will  give 
rise  to  errors,  as  discussed  in  Chap.  II,  Sects.  C  and  D;  and 
Chap «  II I ,  Sect .  D . 

In  the  following  we  apply  the  theory  developed  in  the  pro¬ 
ceeding  sections  to  obtain  the  optimum  waveform' and  receiver 
weighting  function  appropriate  for  mapping  the  lunar  surface  with 
a  radar  of  the  same  carrier  frequency  (440  Me)  and  limit  on 

observation  time  (T  “9  sec),  as  were  reported  for  the 

max 

Millstone  experiment.  The  observation-time  limit  is  set  by  r-f 
oscillator  stability.  From  the  Millstone  data^*  it  is  clear 
that  the  scatterers  on  the  moon's  surface  are  confined  to  a 
range  interval  equal  to  the  moon  radius,  which  corresponds  to  a 
delay  time  or  x  interval  of  23  msec.  Furthermore,  the  scatterers 
are  confined  to  a  range-rate  interval,  at  maximum  libration  rate, 
that  corresponds,  at  440  Me  carrier  frequency,  to  a  Doppler  fre¬ 
quency  or  y  interval  of  about  76  rad  per  sec.  Consequently,  the 
samples  \f/(x  ,2vk/T)  may  be  expeoted  to  be  zero,  or  very  nearly 
zero,  for  x  outside  D  =  23  msec;  and  for  2tt1c/T.  outside  B  -76  rad 
per  sec.  Thus,  a  BD  scatterer  density  is  at  least  approximated; 
and  BD  =  1.7  <  2tt..  Observe  that,  as  long  as  BD  <  2n ,  the  theory 
of  Chap.  II,  Sects.  D,  E,  and  F,  is  applicable;  and  the  period  P 
may  have  any  value  D  <  P  <  2ti/b  or  23  msec  <  P  <  83  msec.  A 
period  of  30  msec  was  used  in  the  Millstone  experiment. 


Suppose  that  a(x),  the  function  that  specifies  the  smoothing 
in  time  delay,  is  chosen  to  be  a  triangular  pulse,  as  shown  in 
rig.  3  where  6  is  the  width  at  the  half-amplitude  points.  The 
function  b(y)  that  specifies  the  smoothing  along  the  Doppler- 
frequenoy  axis  must  be  chosen  so  that  its  Fourier  transform  is 
confined  to  the  interrel  T.  A  suitable  spike-like  form  for  b(y) 
is 

b  (y )  =  si"TfzJ--  (3-65  > 

and  for  this  choice  of  b(y),  the  attainable  y-axis  resolution  is 
limited  to  that  obtaining  for  T  =  T  . 

ilia  X 

Having  chosen  the  smoothing  functions,  we  shall  apply 
directly  the  theory  in  Sects.  E  and  F  to  obtain,  for  any  specific 
P,  6,  and  T,  the  optimum  transmitted  waveform  and  the  corresponds 
ing  receiver  weighting  function.  We  shall  also  investigate  the 
effect  of  changes  in  P,  6,  and  T  on  the  radar  performance. 

From  Fig.  3  and  Eq.  3-55.  we  see  that  the  |  f n | 2  of  the 
optimum  transmitted  waveform  must  follow  a  (sin  x/x)  envelope 
such  that  the  time  autocorrelation  of  f(t)  is  a  periodic  string 
of  triangular  pulses  of  the  form  shown  in  Fig.  3*  Since  the 
phases  of  the  transmitted  harmonics  fp  are  arbitrary,  the  above 
condition  can  be  satisfied  by  a  wide  class  of  transmitted  wave¬ 
forms,  among  which  is  the  periodic  train  of  rectangular  pulses 
of  width  6  and  period  P.  The  relative  simplicity  of  the  trans¬ 
mitter  and  the  receiver  are  important  advantages  of  the  periodic 
pulse  waveform.  However,  other  waveforms,  whose  f  have  differ-, 
ent  phase  angles,  may  result  in  lower  peak  power  requirements 
if  the  conditions  of  the  theory  can  be  met  in  the  construction 
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of  the  transmitter  and  receiver.  The  waveform  actually  used  in 
the  Millstone  experiment  was  an  approximately  rectangular  pulse 
train  with  a  pulse  width  6  =  500  usee  and  a  period  P  =  30  msec. 

It  follows  from  Eqs.  3-59.  3-65.  and  3-64  that  the  appro¬ 
priate  receiver  is  a  matched  filter  with  an  integration  time 
equal  to  T.  For  T  =  =  9  sec,  such  a  receiver  corresponds 

very  closely  to  the  digital  data-processing  scheme  used  in  the 
Millstone  experiment. 

Thus,  if  the  same  smoothing  is  to  be  obtained  as  in  the 
Millstone  experiment,  the  optimum  transmitted  waveform  (aside 
from  the  possibility  of  altering  the  phases  of  the  transmitted 
harmonics  as  discussed  above)  and  the  corresponding  receiver 
weighting  function  obtained  on  the  basis  of  the  dense-scatterer 
theory  are  in  agreement  with  the  waveform  used  in  the  Millstone 
experiment . 

Assuming  that  ^r(x,y)  is  a  rough  scatterer  density,  we  see 

from  Eqs.  3-54,  3“53.  and  3-56  that,  for  a  given  average  power, 

the  average-signal-to-average-noise  ratio  is  independent  of  the 

duration  of  the  observation  interval  T  and  the  smoothing  in 

Doppler  frequency  b(y).  Thus,  Doppler-f requency  resolution  can 

be  increased  without  cost  in  signal-to-noise  ratio  by  increasing 

T  from  any  small  value  to  the  allowable  maximum.  The  interval 

T  =9  sec  used  in  the  Millstone  experiment  may  thus  be  con- 
max 

sidered  optimum  in  view  of  the  equipment  limitations.  In  con¬ 
trast,  as  the  width  of  the  smoothing  function  a(x)  decreases, 
the  average-signal-to-average-noise  ratio  decreases.  The  choice 
of  6  =  500  usee  in  the  Millstone  experiment  is  thus  a  compromise 

between  the  greater  resolution  a  shorter  pulse  would  provide  and 
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the  cost  in  signal-to-noise  ratio  (for  fixed  average  transmitted 
power ) . 

As  long  as  D  <  P  <  2tt/b,  we  note  that:  (l)  the  samples 
^(X,2TTk/T)  are  zero  outside  an  interval  smaller  than  2tt/p  along 
the  y-axis,and  \fl(x,y)  is  zero  outside  an  interval  smaller  than 
P  along  the  x-axis;  and  that  (2)  the  set  of  Gn(x,y)  is  complete 
(assuming  f  f  0  for  all  n)  for  x  inside  an  interval  of  length  P 
and  for  samples  at  y  =  2nk/T  inside  an  interval  of  length  2tt/p. 
Therefore,  if  the  lunar  surface  is  characterizable  by  a  rough 
scatterer  density,  Eq.  3"5^  applies.  As  a  result,  the  average- 
signal-t o-average-noise  ratio  is  independent  of  the  period  P  for 
D  <  P  <  2tt/b.  Furthermore,  the  resolution  is  also  independent 
of  P  for  D  <  P  <  2tt/b. 


CHAPTER  IV 


OOMCU7SZOX8 


A.  GENERAL  THEORY 

There  are  two  basic  limitations  on  radar  performance  in 
dense-soatterer  applications: 

1.  The  incompleteness  of  the  set  of  characteristic  func¬ 
tions  0n(x,y).  This  incompleteness  results  in  an  ambiguous 
component  of  the  scatterer  density,  \f/jjx,y),  which  makes  no 
contribution  to  the  radar  return  and  hence,  even  in  the  absence 
of  noise,  cannot  be  deduced  from  the  received  waveform.  Only 
the  unambiguous  component,  ^||(x,y),  can  be  uniquely  determined 
from  the  noise-free  radar  return. 

2.  In  the  presence  of  additive,  stationary,  white,  Gaussian 
noise,  only  an  estimate  of  a  smoothed  form  of  \f/(x, y)  can  be 
obtained  and,  then,  only  if  the  desired  smoothing  function  is 
expandable  in  terms  of  the  6n(x,y).  As  the  degree  of  smoothing 

is  increased,  the  signal-to-noise  ratio  on  the  radar  display 
Increases;  but  the  sharp  details  of  \f/lx,y)  (or  |\^(x,y)|)  are 
smoothed  out.  Any  radar  design  must  be  a  compromise  between  a 
loss  of  detail  in  excess  of  that  resulting  from  the  incomplete¬ 
ness  of  the  0n(x,y)  and  tbs  signal-to-noise  ratio  on  tbs  display. 
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B.  DELAY-DOPPLER  THEORY 

Time  delay  and  Doppler  frequency  constitute  a  pair  of 
coordinates  for  which  a  scatterer-density  function  can  be  defined 
so  that  its  determination  or  estimation  may  be  taken  as  the 
objective  of  the  radar  system.  f 

If  all  the  time  delays  are  known  a.  priori  to  be  confined  to 
a  finite  interval  D,  and  if  the  observation  time  is  confined  to 
a  finite  interval,  then: 

1.  Only  periodic  waveforms  need  be  considered  as  possible 

transmitted  waveforms. 

2.  A  transmitted  waveform  of  period  P  less  than  D  will 
give  rise  to  an  ambiguous  component  of  the  scatterer 
density . 

3.  Only  the  \j/^(x,y)  component  of  the  scatterer  density, 
that  is,  the  component  that  can  be  expressed  as  a 
sin  x/x  interpolation  of  samples  spaced  at  intervals 
2tt/t  along  y  (see  Eq.  3>ll),  contributes  to  the  radar 
return. 

If  a  periodic  transmitted  waveform  of  period  P  equal  to  D 
is  used,  and  if  it  is  known  a_  priori  that  the  samples  of  »|^p(x, y) 
at  intervals  2tt/t  along  y  are  zero  for  y  outside  an  interval  B, 
then  i//q»(x,y)  can  be  deduoed  unambiguously  from  the  noise~free 
radar  return  if,  and  only  if, 

1.  BD  <  2v 

2.  All  harmonics  of  the  transmitted  waveform  are  different 
from  sero. 

If  the  foregoing  conditions  are  satisfied  and  if  a{x), 
specifying  the  smoothing  along  the  x-axis,  is  a  realizable  auto¬ 
correlation  function,  then,  for  a  given  transmitted  average  power 
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the  output  signal-to-noise  ratio  is  a  maximum  when  the  time 
autocorrelation  funotion  of  the  transmitted  waveform  is  a 
periodic  string  of  pulses  of  the  form  of  a(x).  The  receiver 
weighting  function  appropriate  for  siioh  a  waveform  resembles 
olosely  the  stttehed  filter* 
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Unclassified 
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8/61 

E8L-TM-119 

Unclaiiifted 
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